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Abstract. For the general central force equations of motion in n > 1 dimensions, a com¬ 
plete set of 2n first integrals is derived in an explicit algorithmic way without the use of 
dynamical symmetries or Noether’s theorem. The derivation uses the polar formulation 
of the equations of motion and yields energy, angular momentum, a generalized Laplace- 
Runge-Lenz vector, and a temporal quantity involving the time variable explicitly. A variant 
of the general Laplace-Runge-Lenz vector, which generalizes Hamilton’s eccentricity vector, 
is also obtained. The physical meaning of the general Laplace-Runge-Lenz vector, its vari¬ 
ant, and the temporal quantity are discussed for general central forces. Their properties are 
compared for processing bounded trajectories versus non-precessing bounded trajectories, 
as well as unbounded trajectories, by considering an inverse-square force (Kepler problem) 
and a cubically perturbed inverse-square force (Newtonian revolving orbit problem). 


1. Introduction 

Classical central force dynamics have been long studied from both physical and mathemat¬ 
ical viewpoints. On one hand, many basic physical models are described by central forces, 
e.g. planetary motion and Coulomb scattering, which have inverse-square radial forces; vi¬ 
brations of atoms in a crystal, which have linear radial forces; interactions between a pair of 
neutral atoms or molecules, which have complicated nonlinear radial forces. On the other 
hand, the general equations of motion for central forces have a rich mathematical structure 
and are important examples of completely integrable Hamiltonian systems that possess the 
maximal number of constants of motion in involution. 

The constants of motion for central force dynamics in three spatial dimensions are well- 
known to consist of the energy, the angular momentum vector, and an additional vector 
DEI which is a generalization of the usual Laplace-Runge-Lenz vector known for inverse- 
square central forces [3j. These constants of motion comprise altogether 7 hrst integrals of the 
equations of motion for a general central force. Only 5 of these hrst integrals are independent, 
since the generalized Laplace-Runge-Lenz vector is orthogonal to the angular momentum 
vector, while its magnitude can be expressed in terms of the energy and the magnitude of 
angular momentum, so thus there are two relations among the 7 hrst integrals. The derivation 
of energy and angular momentum can be done in a simple way through Noether’s theorem 
[3] based on symmetries of the Lagrangian formulation of the equations of motion. Noether’s 
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theorem (in its most general form [HIS]) shows that every group of symmetry transformations 
under which the Lagrangian is invariant (to within a total derivative term) gives rise to a 
hrst integral of the equations of motion, and conversely, every first integral arises from a 
group of symmetry transformations of the Lagrangian. Energy is given by the group of 
time-translations, and the components of angular momentum are given by the SO{3) group 
of rotations, where these symmetry groups each act as point transformations on the position 
variables and the time variable. The components of the Laplace-Runge-Lenz vector, in 
contrast, arise from a hidden SO{3) group of symmetry transformations which are not point 
transformations but instead turn out to be dynamical (hrst-order) symmetries when acting 
on the position and time variables. Indeed, a simple, explicit formulation of these hidden 
symmetry group transformations is hard to hnd in the literature |6] and typically only the 
generators are presented BE]. Moreover, while the set of all point symmetries of the central 
force equations of motion is hnite-dimensional, the set of all dynamical symmetries is inhnite- 
dimensional since a general hrst-order symmetry will necessarily involve arbitrary functions 
of all of the hrst integrals [!]• (In particular, for any dynamical system, multiplication of 
any symmetry by any function of hrst integrals automatically yields a symmetry, since hrst 
integrals are constants for all solutions of the system.) Therefore, any direct calculation of 
the hidden SO{3) symmetry group of the central force equations of motion starting from the 
determining equations for dynamical symmetries will require, at some step, the calculation 
of all hrst integrals, which undercuts the entire approach of using Noether’s theorem to 
obtain the Laplace-Runge-Lenz vector. See Ref. [H IH [10] for some alternative approaches 
using symmetry ideas. 

The purpose of the present paper is to take a fresh, comprehensive look at how to derive 
all of the hrst integrals for the general central force equations of motion, without the use of 
dynamical symmetries or Noether’s theorem. This problem will be studied in an arbitrary 
number of dimensions n > 1. Several interesting main new results will be obtained, which 
include deriving a general n-dimensional Laplace-Runge-Lenz vector as a hrst integral in 
two explicit algorithmic ways, and showing a complete set of 2n hrst integrals is generated 
from this general Laplace-Runge-Lenz vector, energy, angular momentum, and a temporal 
quantity involving the time variable explicitly. This temporal quantity provides a hrst in¬ 
tegral that is not a constant of motion. A new variant of the general Laplace-Runge-Lenz 
vector will also be introduced, which arises naturally from these derivations and provides 
a generalization of Hamilton’s eccentricity vector. The physical meaning of the general 
Laplace-Runge-Lenz vector, its variant, and the temporal quantity will be discussed for gen¬ 
eral central forces, and their properties will be compared for processing bounded trajectories 
versus non-precessing bounded trajectories, as well as unbounded trajectories. 

In general, a constant of motion of n-dimensional dynamical equations of motion with 
position variable f and time f is a function C oi r and dr/dt = v whose time derivative 
vanishes 

dC 

^(fW.«W)L.. = o (1.1) 

for all solutions r{t) of the equations of motion. Similarly, a first integral is a function I of 
r, V, and t such that 

^{t,f{t),v{t))L,n,= 
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holds for all solutions r{t). Note that if a first integral does not involve t explicitly, then it is 
a constant of motion. Because dynamical equations of motion in n dimensions comprise a set 
of n second-order differential equations for r{t), the number of functionally-independent first 
integrals is equal to 2n, while the number of functionally-independent constants of motion 
is equal to 2n — 1. 

Our starting point is the well-known fact that [3] every solution r{t) of the central force 
equations of motion in n > 1 dimensions 
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lies in a 2-dimensional plane that is spanned by the initial position vector r(0) and the 
initial velocity vector u(0). Thus the dynamics can be described by polar variables r,6 in 
this time-independent plane. The reduced equations of motion 

— = caV-F — = -2uv/r (1.4) 

comprise a pair of coupled second-order nonlinear differential equations for the dynamical 
variables r(f), 6{t)^ where v = dr/dt is the radial speed and u = d9/dt is the angular speed. 
A hrst integral of the polar equations of motion fll.4|) is a function I{t,r,9,v,u) satisfying 

= It + vIr + ujIe + {uj‘^r + m~^F{r))P-{2uv/r)I^ = 0. (1.5) 

This is the determining equation for all first integrals of these equations of motion. 

In Sec. [21 we will briefly review some aspects of this reduction of the n-dimensional prob¬ 
lem, including the Lagrangian formulation of the polar equations of motion. We also review 
the point symmetries of the polar equations of motion and show that for a general central 
force these symmetries comprise only time-translations and rotations. 

In Sec. 121 we directly derive the general solution of the determining equation fll.5p for 
polar hrst integrals, which yields four functionally-independent hrst integrals. Two of these 
hrst integrals are the energy E and the scalar angular momentum L. The third hrst integral 
is an angular quantity 0 which corresponds to the angular direction of a general Laplace- 
Runge-Lenz vector in the plane of motion. These three hrst integrals E, L, 0 are constants 
of motion. The fourth hrst integral is a temporal quantity T which involves t explicitly. It 
is related to the angular quantity by the property that t = T when 9{t) = 0, for all non¬ 
circular solutions r{t) ^ const.. These two hrst integrals are not widely known or explicitly 
discussed in much of the literature. We conclude Sec. [3] with a detailed discussion of some 
general aspects of the quantities 0 and T, especially their connection to apsides (turning 
points) on trajectories {r{t),9{t)) determined by solutions of the polar equations of motion 
fll.4p . We also discuss a variant of 0 which is connected to other distinguished points on 
trajectories. 

Next, in Sec. 01 we illustrate the physical meaning of 0 and T for two important examples 
of central force dynamics: (1) inverse-square force E = —k/r"^, e.g. Kepler problem; (2) 
cubically perturbed inverse-square force E = —k/r'^ — e.g. Newtonian revolving orbit 
problem (which has precessing trajectories). 

In Sec. O we discuss Noether’s theorem for the polar equations of motion, and we use it 
to derive the symmetries corresponding to the hrst integrals. For the 0 and T hrst integrals, 
these symmetries are found to be dynamical (hrst-order) symmetries which involve both v 
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and cj, while for the L and E first integrals, the symmetries consist of time-translations and 
polar rotations, which are point symmetries. We show that this entire set of symmetries 
generates an abelian algebra (namely, all of the symmetries commnte with each other). We 
derive the corresponding gronps of transformations in an explicit form and work ont how 
each gronp acts on all of the first integrals. In particnlar, these symmetry gronps are shown 
to act in a simple way as point transformations on {t,r,6, L, E). We nse this new resnlt to 
formulate a novel method employing extended point symmetries to derive the first integrals 
0 and T as well as the hidden dynamical symmetry group connected with them. 

In Sec. [6l we express the hrst integrals E, L, 0, T in an n-dimensional geometric form 
in terms of the position vector r and velocity vector n in M”. Through this geometrical 
correspondence, L is shown to yield an antisymmetric angular momentum tensor in the 
plane of motion, and 0 is shown to yield a directional unit vector in the plane of motion. 
We discus the general properties of this unit vector and show how it naturally gives rise to 
a general Laplace-Runge-Lenz vector and a variant vector. 

In Sec. [71 we write out the general Laplace-Runge-Lenz vector and its variant in n dimen¬ 
sions for the two examples of central force dynamics considered previously, and we look at 
the physical and geometrical properties of the resulting n-dimensional vectors. 

Finally, in Sec. [H we make some concluding remarks and outline avenues for future work. 

2. Preliminaries 

The central force equations of motion fll.dp in can be expressed equivalently as a 
hrst-order system 

— =n, — = m^\f\ ^E{\f\)r (2.1) 

for position r and velocity v when n > 1. If these two vectors are collinear for all time t in a 
solution {f{t),v{t)), then the position vector f{t) is clearly conhned to a straight line in M"". 
In this case the motion trivially lies in any time-independent plane that contains this line. If 
instead the vectors r and v are not collinear at some time t in a solution, then span(r, n) at 
time t dehnes a 2-dimensional plane in M”. To see that this plane is time-independent, we 
note that the equations of motion fl2.ll) yield span(dr/dt, dv/dt) = span(F, Ar) = span(F, r^, 
where A = F(|r|)|r|“^. Hence the span is time-independent, which immediately implies that 
the motion lies in a time-independent plane in spanned by the initial position r(0) and 
initial velocity 11(0). 

Let 

{eue2} (2.2) 

be any fixed (time-independent) orthonormal basis for the plane of motion. Then we have 
r = (r cos 9)ei + (r sin 9)e2, v = {v cos 9 — ojr sin 9)ei + {v sin 9 + ujr cos 9)e2 (2.3) 

where r and 9 are polar variables, and v = dr/dt and u = d9/dt are their time derivatives. 
Substitution of this polar representation into the equations of motion fl2.ip yields the polar 
equations of motion fll.dp in which r and 9 are the dynamical variables. This reduces the 
number of degrees of freedom from n > 1 to 2. 

The polar equations of motion fll.4p arise from the well-known Lagrangian 

L = |m(n^ -f r^w^) — U{r) 
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where U{r) is the radial potential determined (up to an additive constant) by the central 
force 


F(r) = -U'{r). 


(2.5) 


In particular, we have 


6L 

6r 


mu'^r + F{r) — m 


(Pr 

dP ’ 


5h 2^20 

— = —2muvr — mr 

59 dP 


( 2 . 6 ) 


Hereafter, we will put m = 1 without loss of generality (via rescaling the physical units of 
the dynamical variables). 

In polar variables, a point transformation is an invertible mapping of the coordinate space 
{t, r, 6) into itself. We will be interested in Lie groups of point transformations. A one¬ 
dimensional Lie group of point transformations, with group parameter e, consists of 


f ^ tp)(t, r, 0), r —)■ r(e)(f, r, 6^), 9 ^ 9(^){t,r,9), oo < e < oo 


(2.7) 


such that e = 0 yields the identity transformation. The generator of the group is an inhni- 
tesimal transformation 


5t = 


dt 


h) 


de 


e=0 


= r(t, r, 9), 5r = 


dr 


h) 


de 


e=0 


= ^{t,r,9), 59 = 


d9 


id 


de 


e=0 


= P{t,r,9) (2.8) 


which can be viewed as dehning a vector held 


X = Tdt + idr + pde (2.9) 

on the space of variables (f, r, 9). From this inhnitesimal generator, the point transformations 
(12.7p can be obtained by exponentiation of the vector held (12.9p . 

Point transformations have a natural action on functions (r(t), 9{t)). In inhnitesimal form, 
this action is given by 


ir{t) = {((, r{t),e(t)) - v{t)T{t, r(f),«(()), Se{t) = »(f)) - u(()t((, r(f),«(()) 

( 2 . 10 ) 

where v{t) = dr(t)/dt and oj(t) = d9(t)/dt. The corresponding generator can be expressed 
as a vector held 

± = P^dr + P^dg, P^ = ^-VT, P^ = p-UJT. (2.11) 

This vector held is called the characteristic form m associated to the generator fl2.9p . 

Both vector helds fl2.1ip and fl2.9p can be prolonged to act on time-derivatives of the polar 
variables. The prolongation of the vector held fl2.1ip is very simple, 

dP'’ dP® 

prX = X H—H—+ • • • (2-12) 

dt dt 

while the prolongation of the vector held (12.9p takes the related form 

prX = prX-I-r-f- (2.13) 

dt 

where the total time derivative is viewed as a vector held 

— = dt P vdj. P ujdg (2.14) 

dt 

by the chain rule. Note these prolonged vector helds are dehned in the coordinate space 
{t, r, 9, v,u,...), which is called the jet space of the polar equations of motion. 
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A point symmetry of the polar equations of motion fll.4p is a Lie group of point transforma¬ 
tions fl2.7D that leaves invariant the solution space of the equations. Inhnitesimal invariance 
of the solution space is expressed in terms of the prolonged generator fl2.9p of the transfor¬ 
mation acting on the polar equations of motion by 


prX 


(fr 

dt'^ 


— u'^r — F{r) 


soln. 


= 0, prX 


d‘^9 

dF 


+ 2uvlr 


= 0 


soln. 


(2.15) 


which is required to hold for all solutions of the equations. This invariance condition can be 
expressed more simply through the characteristic generator fl 2 . 1 ip . which has an equivalent 
action on the equations of motion because of the relation fl2.14p . In particular, since a total 
time-derivative of the equations of motion necessarily vanishes on solutions of the equations, 
the invariance condition becomes 


d^e 


■ uP'r 


F{r) 


soln. 


d^P^ 

dF 




dP^ 

2ojr -F'(r)P" = 0 

dt 


prX ( — 4 - 2uv/r 


soln. 


d^p^ dP^ 
+ 2 


dF 


dt 


dP^ 

vr~^ + 2ur~^—r. - 2ojvr~‘^P^ = 0 


(2.16) 


dt 


where d?r/dF, d?9/dF, and all higher-order time derivatives are eliminated through use of 
the equations of motion fll.4p . Thus, the generator of a point transformation fl2.9p will be 
an inhnitesimal symmetry of the polar equations of motion iff the associated characteristic 
generator fl2.1ip satishes the invariance condition fl2.16p . For this reason the two generators 
fl2.9p and fl 2 . 1 ip are commonly referred to as equivalent symmetry vector helds, and the 
invariance condition fl2.16p is called the determining equation [H |5] for point symmetries. 

The determining equation fl2.16p is a straightforward linear partial differential equation 
to be solved for the inhnitesimal symmetry components r(t,r, 6 *), ^(t,r,9), and 
In particular, this equation will split with respect to v and u in the jet space (t,r,9,v,oj), 
yielding a linear overdetermined system on the functions T(t,r,9), ^(t,r,9), and 6 *). 

If we regard the force F{r) as an additional unknown, then we can also determine any 
inhnitesimal point symmetries that hold only for special force expressions F{r). In all cases, 
the set of admitted inhnitesimal point symmetries forms a Lie algebra. 

By a direct calculation (using Maple), we obtain the following result. 

Theorem 1. The infinitesimal point symmetries 02.91) admitted by the polar equations of 
motion 01.4p with a nonlinear force F{r) are spanned by the generators: 

(a) general F{r) 

Xi=a 0 X 2 = ai, 

(b) F{r) = kr^, p 7 ^ 1 

Xi = de, X 2 = dt, X 3 = tdt + -^rdr 

1 — p 

(c) F{r) = kr + kr~^, k 0,k 0 

Xi - de, X 2 — dt, X 4 — exp{2k^:H) (dt + k^^^rdr) 


(2.17) 

(2.18) 
(2.19) 
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The point symmetry transformations generated by Xi, X 2 , X 3 , and X 4 are, respectively, 
given by 


( 1 ) 9 ^ 9 + e polar rotation 

(2) t —)• t + e time-translation 

(3) t —)■ exp(e)t, r —)■ exp (2e/(l — p)) r scaling 


( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 


(4) t ^t + k-^/Hn{{l-2k^/^eexp{2k^/H)) r ^ r {l - 2k^/^eexp{2k^/H)) 


(2.23) 


time-dependent dilation 


Notice that only two point symmetries are admitted for a general central force F{r), 
whereas the equations of motion possess four functionally independent hrst integrals. Hence, 
point symmetries of the equations of motion are not rich enough to capture all of the hrst 
integrals. 

An important generalization of point symmetries is provided by dynamical symmetries. A 
dynamical symmetry |T2] of the polar equations of motion fll.dp consists of an inhnitesimal 
transformation of the general form fl2.9p in which the components r, ip are allowed to 
depend on the radial speed v and angular speed u (in addition to the variables t, r, 9) such 
that the invariance condition fl2.15p holds for all solutions of the equations of motion. As 
in the case of point symmetries, it is simpler to work with the equivalent characteristic 
generator fl2.1ip for dynamical symmetries, which satishes the determining equation fl2.16p . 

Similarly to point symmetries, dynamical symmetries can be exponentiated to obtain 
a group of transformations acting on solutions of the polar equations of motion. These 
transformations take the form of point transformations acting on the dynamical variables 
{r(t),9(t),v(t),oj(t)) given by any solution {r(t),9(t)), but they cannot be extended off so¬ 
lutions to act on the coordinate space (t,r,9,v,oj). 

Moreover, in further contrast to point symmetries, the determining equation 02.161) for 
dynamical symmetries no longer splits with respect to v and u in the jet space (t, r, 9, v, co), 
since these variables now appear in the dynamical symmetry components T{t,r,9,v,uj), 
f(t,r,9,v,oj), and 'ip{t,r,9,v,oj). This means that the determining equation is a coupled 
pair of linear partial differential equations, and in general we cannot solve these equations 
without already knowing how to integrate the polar equations of motion themselves. In par¬ 
ticular, the general solution of the dynamical symmetry determining equation will involve 
arbitrary functions of all hrst integrals of the polar equations of motion. 

Therefore, it will be more useful to derive the hrst integrals directly from the equations 
of motion, rather than use a symmetry approach. Nevertheless, dynamical symmetries have 
a direct connection to hrst integrals through Noether’s theorem, which can be used to hnd 
the specihc dynamical symmetries that correspond to the hrst integrals in Theorem [2J 


3. Derivation and properties of polar first integrals 


For the polar equations of motion fll.4l) of central force dynamics, all hrst integrals are 
functions I(t,r,9,v,u) determined by equation fll.Sp . With m = 1 and F(r) = —U'{r), this 
determining equation becomes 


dl 

dt 


It + vir ule + (ca^r - U'{r))I„ - {2uvr = 0 


soln. 
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which is a linear first-order partial differential eqnation. Its general solution can be obtained 
through the method of characteristics HU by integrating the system of differential equations 

dt dr dO dv du dl 

1 V 


V u u‘^r — U'{r) —2uvr~^ 

This system can be arranged in a triangular form 

du —2u 
dr r 
dv u'^r — U'{r) 

dr V 

dd u 

dr V 
dt 1 

dr V 
dl 


0 


dr 


= 0 


(3.2) 


(3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 


whereby the successive integration of these differential equations fl3.3l) - (l3.6p is achieved by 
separation of variables. Each constant of integration will then be a particular first integral 
satisfying the determining equation fl3.ip . 

Separating variables in equation fl3.3p . we obtain the first integral 

Ji = (3.8) 


We can now express u in terms of Ji and r: 

u = 

Then equation fl3.4p is separable 

dv Ii^r~^ — U'{r) 
dr V 

Hence we obtain a second first integral 

/2 = l^;2+ l/2^-2 + 17(r). 

We can next express the magnitude of v in terms of Ji, J 2 , and r: 

\v\ = V 2(/2 - u{r)) - h^r-^. 

Both equations (13.5p and (13. 6 p now become separable, 

dO sgn(n)Ji 


dr r2^2(J2 - U{r)) - h^r-'^ 


and 


dt 


sgn(n) 


dr ^2{h-U{r))-h^r-^' 
Hence we obtain two more first integrals 

f sgn(n) 


v/2(/2 - U{r))r^ - h^r^ 


dr = Q 


(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 


h = e-h 





















and 


L = t — 


sgn(n) 


VW2 


■.dr = T. 


(3.16) 


f/(r)) - /i 2 r -2 

Finally, the remaining differential eqnation fl3.7p shows that I is an arbitrary fnnction of 
the previons fonr first integrals fl3.8p . fl3.1ip . fl3.15p . fl3.16p . The following resnlt is now 
immediate. 


Proposition 1. The general solution of the determining eguation fl3.ip for first integrals of 
the polar eguations of motion fll.4l) is given by I = /(/i, J 2 , Is, 14 ) where f is an arbitrary 
differentiable function and where Ii, I 2 , I 3 , Ia have physical units of angular momentum, 
energy, radians, and time, respectively. 

Note that these first integrals Ji, I 2 , h, I a are fnnctionally independent becanse they each 
have different physical nnits. Hence, Ji, J 2 , I 3 , I a provide the complete qnadratnre of the 
polar eqnations of motion. 


3.1. Normalization (“zero-point” values). Any first integral I remains conserved if an 
arbitrary constant is added to it. This freedom represents the choice of a “zero-point” valne 
for the physical qnantity defined by the first integral. In particnlar, we can write 


Ji — L Lq (3.17) 

I 2 = E -\- Eq (3.18) 

/3 = 0 + 00 (3.19) 

/4 = T + To (3.20) 


where Tq, Eq, 0o, Tq will denote the zero-point constants, and T, E, 0, T will be the 
normalized physical qnantities, which have nnits of angular momentum, energy, radians, 
and time, respectively. 

The zero-point constants in the first integrals need to be specified by some additional 
considerations such that the resulting quantities L, E, 0, T are physically meaningful and 
mathematically well-defined for all solutions of the polar equations of motion fll.4p . It is clear 
how we can do this for the angular momentum constant Lq and the energy constant Eq, by 
adopting the standard Newtonian expressions for the physical angular momentum L and the 
physical energy E for motion under a central force. But since the physical interpretation of 
0 and T is not obvious just from their expressions, we will instead use a different argument 
to determine all four constants Lq, Eq, 0o, Tq based on general properties of the effective 
potential for the equations of motion. Normalization of 0 is often overlooked in the literature 
yet is crucial for understanding its relationship with the Laplace-Runge-Lenz vector. 

From the first integrals fl3.8l) and (13.111) . written in the respective forms fl3.17|) and fl3.18p . 
the effective potential is defined by 

17eff.(r) = \{L + To)V-2 + U{r) - Eq. (3.21) 

This potential determines the types of trajectories admitted for solutions of the polar equa¬ 
tions of motion (11.4p . In all cases of physical interest, we may suppose that the effective 
potential has at least one eguilibrium point, r = req., defined by the condition that the central 
force —U'{r) vanishes at r = req, (which can include r = 0 or r = cx)). This condition coin¬ 
cides with the effective force —Heff (^) being zero when both the radial and angular speeds 
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are zero. Thus, the set of equilibrium points is given by the roots of the equation 

0 ~ (’"eq.) = (^eq.) (3.22) 

The physical meaning of an equilibrium point is that it corresponds to a static solution of 
the polar equations of motion fll.dp . with r = r^q. = const, (and 6 = const.). 

A natural condition to determine Lq and Eq is that both L\j.=r^^ and E\r=r^^ must be 
zero when these two hrst integrals are evaluated for static solutions of the polar equations of 
motion fll.dp given by an equilibrium point fl3.22p away from the origin, r = req. 7 ^ 0. Since 
a static solution has a; = n = 0 , the condition L\r=r^^^ = 0 directly implies 

To = 0 (3.23) 

in the hrst integral fl3.17p . while the condition E\r=req = 0 together with L\r=req = Lq = 0 
yields 

Eo = U{req.) (3.24) 

in the hrst integral fl3.18p . This gives 

L = ur^ (3.25) 

and 

E=y + \L^r-^ + U{r) - U{r,q) (3.26) 

which are the usual Newtonian dehnitions of physical angular momentum and physical energy 
for central force dynamics. Note the ehective potential now simplihes: 

17eff.(r) = iLV-2 + U{r) - U{r,q). (3.27) 

This equilibrium point argument, however, will not extend directly to ©o and Tq, because 
0 and T can have any values at an equilibrium point r = rgq,. To determine ©o and Tq, 
we will consider, more generally, distinguished points of the ehective potential other than 
equilibrium points. Two natural distinguished points are inertial points and turning points, 
as these points are dehned solely in terms of E. 

An inertial point is a hnite radial value r = r* at which the ehective force vanishes. The 
set of all inertial points is thus given by the roots (if any) of the ehective force equation 

0 =(r*) = LV-3 - f/'(r*), 0<r*<cx). (3.28) 

This is a generalization of the equation fl3.22p dehning an equilibrium point ifL 7 ^ 0. Hence, 
a solution of the polar equations of motion fll.41) possesses an inertial point if (and only if) 
both dv/dt = 0 and a; = L/r^ 7 ^ 0 hold at a point on the trajectory where r = r*. 

A turning point is a hnite radial value r = r* at which the ehective potential is equal to 
the energy E. Hence the set of all turning points is given by the roots (if any) of the energy 
equation 

0 = Heff.(n) - 7? = + t/(r*) - E - H(req.), 0 < r* < 00 . (3.29) 

A solution of the polar equations of motion (ll.4p possesses a turning point if (and only if) 
n = 0 holds at a point on the trajectory where r = r*, since E — UesX^) — from equations 
fl3.26p and fl3.27p . Note that a turning point r = r^ will coincide with an inertial point r = r* 
whenever the energy has the value E = Hefr, (’"*)• 

For all potentials U{r) of physical interest, at least one turning point or one inertial point 
can be assumed to exist for every solution of the polar equations of motion (II.4p . However, if 
a solution is circular, then the integral expressions (I3.19p and (I3.20p cannot even be dehned 
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since r is constant, and therefore only the hrst integrals for angular momentum fl3.25p and 
energy fl3.26p exist. Consequently, the question of determining the constants ©o and Tq 
needs to be considered only for non-circular solutions. 

A natural way to determine ©o and Tq is by requiring that the respective values of the hrst 
integrals fl3.19p and fl3.2Up for each non-circular solution of the polar equations of motion 
(El are given by the conditions ©|r-=ro = and T\r=ro = to where either r(to) = tq = r* 
and 0(to) = Oq is an inertial point on the trajectory, or r(to) = tq = r* and 6*(to) = 6*o is a 
turning point on the trajectory. These conditions lead directly to the expressions 


and 


where 


Q = e-L 


sgn(t;) 


T = t- 


Iro ^2{E + U{r,J-U{r))r^-L^r^ 
r sgn(n) 


dr 


'ro 


^2{E + 17(req.) - U{r)) - L^r-^ 


dr 


ro = r^, rl{U{r^)-E-U{r^^)) = \L^ 


(3.30) 

(3.31) 

(3.32) 


or 

ro = r*, r*^U'{r*) = L\ (3.33) 

Note that the choice of rg in these integrals corresponds to specifying the values for the 
constants ©o and Tq. 

Thus we have proved the following main result. 


Theorem 2. For the polar equations of motion (II.4p of general central force dynamics: 

(1) L and E are well-defined first integrals for all solutions. (L depends solely on r,uj; E 
depends solely on v and Ues.{T)). 

(2) © and T are well-defined first integrals for all non-circular solutions. (Both © and T 
depend solely on L, E, sgn(n) and Uesfr).) 

(3) L, E, © are functionally independent constants of motion. 

(4) Every first integral is a function of L, E, ©, T, and every constant of motion is a function 
ofL, E, ©. 


We emphasize that these four hrst integrals, given by equations f|3.25|) . fl3.26p . fl3.30p . 
fl3.3ip . can be directly verihed to obey 

dL I dE I d© I dT I ^ 

We will see later that the hrst integral © is directly related to the angle of the Laplace- 

Runge-Lenz vector in the plane of motion for non-circular solutions. 


3.2. Evaluation using turning points. Each solution {r{t),6{t)) of the polar equations 
of motion fll.4p describes a trajectory in the plane of motion. The shape of a trajectory 
is a curve 6 = /(r, L, E, ©) parameterized by the values of the angular momentum L and 
the energy E, as well as the angle ©, with r belonging to some specihed radial domain 

^min Fi Fi ^max- 

The angular hrst integral fl3.30p gives an algebraic (quadrature) equation for the curve. 


e = f{r,L,E,Q) = Q + L 


sgn(n) 


Iro V2(E + 17(req.)-f/(r))r4-LV 
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dr. 


(3.35) 






























while the motion along the curve is implicitly given by the temporal first integral fl3.3ip . 


t = T + 


sgn(u) 


'ro 


^2{E + U{r,J-U{r))-L^r-^ 


dr. 


(3.36) 


If the angular momentum L = r'^u is equal to zero, then the trajectory is purely radial, 
since cu = 0 implies 9(t) = 9(0) is constant. The equation of the curve (13.351) thereby reduces 
to an unbounded radial line 9 = 0 = const, through r = 0, where the angular first integral 
fl3.30p evaluated on the radial trajectory (r(t),9(t)) is given by 0 = 6^(0), independently of 
the choice of tq. Thus, in this case, 0 is the angle of the radial trajectory in the plane of 
motion. 

The physically more interesting trajectories are non-radial, which occur if the angular 
momentum has a non-zero value, L = r^cj ^ 0, so then 9(t) is no longer constant. There are 
two different primary types of non-radial trajectories: bounded and unbounded. A bounded 
trajectory has a radial domain given by 0 < rmin < rmax < oo, whereas the radial domain 
for an unbounded trajectory is given by 0 < rmin < r^ax = oo. A crucial feature of both 
types of trajectories is the number of turning points that occur in the radial domain. When 
a trajectory has rmin > 0 , any point with r = r^^^i on the trajectory corresponds to a 
turning point given by r* = rmin. Similarly, when a trajectory has r^ax < oo, any point with 
r = Tmax on the trajectory corresponds to a turning point given by r* = rmax- Bounded non- 
radial trajectories are thus characterized by a radial domain having precisely two turning 
points r* = and r* = Tmax when the trajectory is non-circular, or a single turning point 

= 'Train = ’"max in the Special case when the trajectory is circular. Unbounded non-radial 
trajectories are characterized by having a radial domain with either only one turning point, 
T* = Tmin, or no tuming points. 

Any turning point that occurs on a non-circular trajectory (r(t), 9(t)) is a radial extremum 
which is either a periapsis of the trajectory where r = r^, = rmin > 0, or an apoapsis of the 
trajectory where r = r^ = r^ax < oo. These apses may occur at any number of distinct 
angles (with the same radial distance r = r*) depending on the shape of the trajectory. When 
a non-circular trajectory (r(t),9(t)) passes through any apsis, (r*,^*), the radial speed v(t) 
will change its sign. From equation fl3.35p we then have the following result. 


Proposition 2. For any solution of the polar equations of motion (11.4p yielding a non¬ 
circular trajectory, the curve (I3.35P determined by the shape of the trajectory is locally sym¬ 
metric (under reflection) around the radial line connecting the origin r = 0 to any apsis 
point in the plane of motion. 


We will refer to the radial line 9 = 9^ determined by a given apsis (r*, 6 **) on a non-circular 
trajectory as an apsis line. Note that every apsis on a given trajectory is a turning point of 
the effective potential, but the set of turning points (I3.29p may include points that do not 
occur on the trajectory. 

It is clear that turning points are important in the evaluation of the angular and temporal 
first integrals (l3.3Up and (I3.3ip since the expression 

^J2{E + U(req,) - U(r)) - L'^r-'^ = |u| (3.37) 

appearing in these integrals will vanish at all turning points r = r^ that occur on a given 
non-circular trajectory (r(t),9(t)). 


12 





















3.3. Piecewise property (“multiplicities”) and trajectory shapes. There are three 
different types of non-circular trajectories which may arise, depending on the number of apis 
points. The physical and mathematical properties of the angular and temporal first integrals 
fl3.30p and fl3.3ip differ in each case. 

First, suppose a non-circular trajectory {r(t),9(t)) possesses no apsis. This implies that 
the radial domain of the trajectory contains no turning points and hence the trajectory 
is unbounded. Consequently, the equation fl3.35p of the curve describing the shape of the 
trajectory will be valid on the whole radial domain, and similarly the angular hrst integral 
fl3.30p will give a unique value for 0 when it is evaluated on any part of the trajectory, 
regardless of the choice of rg in the integral. The most physically meaningful and mathe¬ 
matically simple choice for rg will be an inertial point, tq = r*, given by either the maximum 
or the minimum extremum of the effective potential. With this choice, the value of 0 will 
be the angle 6 = 6* = Q of the point on the trajectory at which the radial speed v is an 
extremum, while the value of T given by the temporal hrst integral fl3.3ip will be the time 
t = t* = T at which this point is reached on the trajectory. 

Next, suppose a non-circular trajectory possesses a single apsis point, {r,9) = (r*, 6 '*). 
The radial domain therefore contains exactly one turning point, r = r^, and the trajectory is 
thus unbounded. In this case the hrst integrals fl3.30l) and fl3.3ip are most naturally dehned 
by having ro = r* chosen to be the turning point (corresponding to the apsis). Then the 
equation fl3.35p of the curve which describes the shape of the trajectory will be valid on 
the whole radial domain. Because this domain contains a single turning point, r = r*, the 
curve will be globally rehection-symmetric around the corresponding apsis line 9 = 9^,, and 
the resulting two halves of the curve will come from separately putting sgn(u) = -|-1 and 
sgn(u) = — 1 in equation fl3.35p . Moreover, the angular hrst integral fl3.3Up will give a unique 
value for 0 when it is evaluated on any part of the trajectory. The physical meaning of this 
value is that it will be the angle 9 = 9^ = Q at which the apsis occurs on the trajectory in 
the plane of motion. The same argument shows that the temporal hrst integral fl3.3ip gives 
a unique value for T which will be the time t = T at which the apsis r = r* is reached 
on the trajectory. 

Finally, suppose a non-circular trajectory {r{t),9{t)) possesses multiple apsis points (ap¬ 
sides). The radial domain of the trajectory thereby contains exactly two turning points, 
which are the endpoints of the domain. Correspondingly, the trajectory is bounded. The 
equation fl3.35p of the curve then will divide up into separate pieces that are determined in 
the following way by the apis points. Since r{t) is increasing nearby any periapsis point, the 
trajectory must pass through an apoapsis point before reaching another periapsis point. Like¬ 
wise, since r{t) is decreasing nearby any apoapsis point, the trajectory must pass through 
an periapsis point before reaching another apoapsis point. Hence, multiple apsis points 
come in pairs consisting of both a periapsis and an apoapsis. Let (r, 6*) = (r*!,^*!) and 
(r, 9) = (r* 2 , 6 ** 2 ) be any pair of adjacent periapsis and apoapsis points on the trajectory, 
whereby 0 < rmin = r*i < r *2 = Tmax < 00 . Consider the piece of the curve starting at the 
periapsis point, as dehned by choosing tq = r*i in the curve equation fl3.35p . with 0 = 0 * 1 . 
This piece of the curve has the radial domain r*i < r < r* 2 , where the other endpoint r = r *2 
is an apoapsis point. Since the curve is locally rehection-symmetric around the apoapsis line 
9 = 0*1 in the plane of motion, the next piece of the curve is dehned by choosing tq = r *2 
in equation fl3.35p . with 0 = 0 * 2 , where the radial domain is r *2 > r > r*i. Note that 
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sgn(t>) differs on these two pieces of the curve. Clearly, this process can be continued to 
piece together the entire curve, which yields the complete trajectory {r(t),9(t)). 

Thus, for any bounded non-circular trajectory with multiple apsis points, it follows that 
the angular hrst integral 03.301) will be multi-valued when it is evaluated on different parts 
of the trajectory, corresponding to the piecewise composition of the curve 03.351) . The same 
property will hold for the temporal hrst integral 03.311) . since it is uses the same value of tq 
that is chosen in the angular the hrst integral 03.3Up . Each pair of values 0 and T given by 
these hrst integrals are the angle and the time at which the trajectory reaches each apsis. 

As a consequence of Proposition [2], the angular separation between two successive apsis 
points on a bounded non-circular trajectory is the same on all parts of the trajectory. From 
the curve equation 03.35|) . it is straightforward to get an integral expression for this angular 
separation. 


Proposition 3. For any solution of the polar equations of motion fll.4p yielding a bounded 
non-circular trajectory, the angular separation between any two successive periapsis points 
or apoapsis points on the trajectory is given by 


Ae = 2L / , 

V2(i? + f/(req.)-17(r))r4-LV2 

and the corresponding time interval is given by 

At = 9 r"“_ SEn(^) 

A,. y2(B + - U{r)) - LV-2 

These expressions are related to the first integrals 0 and T by 

= sgn(t;)(0|^,=^^^^ - 0|r-o=r,„iJ, = sgn(n)(T|^„=^^^^ - 


dr {modulo 27i), 


dr. 


(3.38) 


(3.39) 


(3.40) 


If the angular separation A6 is a rational multiple of 27r, then the angular hrst integral 
fl3.30p will yield a hnite number of distinct values 0 modulo 27r. In this case, the curve 
describing the shape of the trajectory is closed. Note that the apsis points of the trajectory 
in the plane of motion are precessing unless A6 is exactly equal to 27r. 

In contrast, if the angular separation A6 is an irrational multiple of 27r, then the angular 
hrst integral fl3.30p will yield an inhnite number of distinct values 0 modulo 27r. This means 
that the curve describing the shape of the trajectory is open, such that the apsis points of 
the trajectory are precessing in the plane of motion. 

For both kinds of (bounded non-circular) trajectories, the temporal hrst integral 03.311) 
yields a periodic inhnite sequence of values T. 


4. Examples 

The results in Theorem [21 Proposition [2] and Proposition [31 will now be illustrated for two 
examples of central forces. 

4.1. Inverse-square force. Consider the central force 

F = -kr-"^ (4.1) 


U{r) = —kr ^ 
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with the potential 


(4.2) 
























UeS.{r) = \L‘^r - kr L 7 ^ 0 . 

(4.3) 

E = ^v‘^ + \L‘^r-^ - kr-^ 

(4.4) 


This force will be attractive if fc > 0, e.g. planetary motion, or repulsive if A; < 0 , e.g. 
Coulomb scattering of charged particles. In either case, the potential has only one equilibrium 
point, Teq. = oo. Siuce fA(req.) = 0, the effective potential is given by 


Thus 


is the energy first integral fl3.26p . All trajectories can be classified qualitatively from the 
equation E = 

In the repulsive case fc < 0, has no extremum, and E is non-negative. Hence all 

trajectories are unbounded. 

In the attractive case fc > 0, Ues.{'^) has one extremum, which is a negative minimum, 
= —\{k/L)‘^. Hence, for L 7 ^ 0, all trajectories with E > Q are unbounded, while all 
trajectories with 0 > are bounded, and bounded trajectories with E = are 

circular. 

We now evaluate the angular and temporal first integrals fl3.30p and fl3.3ip for the case 
/c > 0, L 7 ^ 0, i.e. the Kepler problem. We will need the relation 

|n| = \/2E — + 2kr~^. (4.5) 

4.1.1. Turning points for the Kepler problem. The turning points r = r* of the effective 
potential fl4.3p are obtained from the energy equation fl3.29p . This is a quadratic equation 


0 = 2Er‘^ + 2kr — 


(4.6) 


with the discriminant D = 4{k^ — 2EL‘^). Turning points exist only when H > 0, where 
D = 0 determines the minimum energy 

k^ 


E^mn 


2L2 


< 0 . 


(4.7) 


Since the trajectories with E = i^min are circular, only the trajectories with E > T^min need 
to be considered (as the first integrals 0 and T exist only for non-circular trajectories). 

For 0 > E > the turning points are given by 

fc ± - 2|E|L2 


K± = 


2\E\ 


(4.8) 


The angular hrst integral fl3.30l) with uq = can be straightforwardly evaluated to yield 


Q± = 0 -h arctan 


— kr 
Lvr 


± sgn(r>L) 


TT 


(4,9) 


after use of the relation (14.5p . Similarly, the temporal first integral (I3.3ip with ro = 
yields 


T+ = t 


rv 


+ 


k 


arctan 


/ 2\E\r — k 

[ VW\v 


TT 


Tsgn(v)-). 


(4.10) 


2\E\ ' (2|E|)3/2 

From the first integral fl4.9p . trajectories with > E > F^min are algebraically described by 

r{e) = 


k T ^/k^-2\E\L^ cos( 0 ± - 9) 
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curves 


(4.11) 



























on the angular domain 


(4.12) 


- + Q±<e < + 0 ± 

where, as shown by Proposition [3l 

= sgn(n)(0_,_ — 0_) = TT (4.13) 

is the angular separation between the two turning points on the curve. These curves are 
ellipses [3], which are bounded and have r^ax = 'f*+ and = r*_, such that r = 0 is one 
focus of the ellipse in the plane of motion. Moreover, in accordance with Proposition |21 each 
ellipse is reflection-symmetric around a radial line connecting the origin r = 0 to either of 
the points r = r*±. The resulting two pieces of the ellipse are each dehned on the radial 
domain rmm <r< r^ax and differ by sgn(v) ^ 0 (where sgn(n) = 0 at each turning point). 

On each elliptic trajectory, both hrst integrals 0± and T± are piecewise functions. Their 
physical and geometrical meaning depends on the choice of the radial value tq = r+ or 
To = r_. If tq = r+ then (r(T+), 6{T^)) = (r^ na x, 0+) is an apoapsis point on the trajectory, 
or instead if tq = r_ then {r(T_),9{T_)) = (rmm, 0-) is a periapsis point. In both cases, 
Q± and T± are continuous at the apsis point r = r*±, but have a jump discontinuity at the 
opposite apsis r = The jump in 0± is equal to 

= 27r (4.14) 


which is the angular separation fl3.38|) between two successive periapsis points or two suc¬ 
cessive apoapsis points on the trajectory. Since A^ is exactly 27r, these trajectories are not 
processing, and hence both of the angular hrst integrals 0 ± are single-valued modulo 27r. 
The jump in T± is simply the period of a single closed orbit. 


At = 


nk 




(4.15) 


Both of the temporal hrst integrals T± thus yield a sequence of times such that the time 
interval (13.391) is the period At. 

Therefore, when the angular hrst integral 0_ and temporal hrst integral T_ are evaluated 
at any point {r{t),9{t)) on an elliptic trajectory (14.lip , we obtain the angle of the radial 
line that intersects the periapsis point and the time at which this point is reached on the 
trajectory. Likewise the angular hrst integral 0+ and temporal hrst integral T+ evaluated 
at any point {r{t),9(t)) yield the angle of the radial line that intersects apoapsis point and 
the time at which this point is reached on the elliptic trajectory. 

Next, for E = 0, there is just a single turning point 


r* = 


2k' 


The angular hrst integral (I3.30p with ro = r* is straightforward to evaluate, 

' LF' — kr" 


Q = 9 + arctan 


Lvr 


sgn(r;L) 


TT 


(4.16) 


(4.17) 


Similarly, evaluation of the temporal hrst integral (13.3 Ih with ro = r* yields 

rv{L‘^ — kr) 


T = t- 


3P 
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(4.18) 
















From the first integral fl4.17p . trajectories with = 0 are algebraically described by curves 


r{9) = 


L 2 


k{l + cos(0 — 9)) 


(4.19) 


on the angular domain 

-7r + 0<0<7r + 0. (4.20) 

These curves (I4.19p are parabolas [3] with the focus at r = 0, such that rmin = r*. In 
accordance with Proposition [2l each parabola is reflection-symmetric around the radial line 
connecting the origin r = 0 to the periapsis point r = r* on the trajectory. The resulting 
two pieces of the parabola are each dehned on the radial domain Tmin < r < oo and differ by 
sgn(n) ^ 0 (where sgn(r>) = 0 at the turning point). 

The hrst integrals 0 and T are piecewise functions on each parabola, such that 
{r{T),9{T)) = (^1^111,0) is the periapsis point. Both 0 and T are continuous at this point. 
Thus, when these hrst integrals are evaluated at any point (r(f), 9{t)) on a parabolic trajec¬ 
tory 04.191) . we obtain the angle of the radial line that intersects the periapsis point and the 
time at which this point is reached on the trajectory. 

Last, for > 0, there is again a single turning point 


r* = 


^^WT2EI?-k 
2E ■ 


(4.21) 


Evaluation of the angular hrst integral O3.30p with tq = r* straightforwardly yields the same 
expression (I4.17p as in the case E = 0. The temporal hrst integral 03.3ip with tq = r* gives 


k 


T = t - (rv + , _ 

2E V 


arctanh 


\/2E' 


vr 


2Er + k 


(4.22) 


From the hrst integral 04.17P combined with the relation 04.Sp . trajectories with E > 0 are 
algebraically described by curves 

(4.23) 


k + yW+2EI? cos{Q — 9) 

— arctan (^2E\L\/1^ + Q < 9 < — + arctan (^2E\L\/1^ + 0. 


(4.24) 


on the angular domain 

TT 

~ 2 

These trajectories are hyperbolas [3], which are unbounded and have r = 0 as the focus. 
The properties of the angular hrst integral 0 and temporal hrst integral T for a hyperbolic 
trajectory are the same as in the parabolic case. 

4.1.2. Inertial points for the Kepler problem. The inertial points r = r* of the ehective 
potential 04.3p are obtained from the ehective force equation 03.281) . This reduces to a linear 
equation 

0 = - fcr (4.25) 

giving 




r 
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(4.26) 





























From Proposition [2l every trajectory with L 7 ^ 0 will contain two inertial points r = r*, 
differing by sgn(r;) ^ 0. In particular, the radial speed at the these points is given by 


V = 




(4.27) 


by combining equations (I4.26p and (14.hh . 

For 0 > > F^min, the angular hrst integral (13.301) with tq = r* can be straightforwardly 

evaluated to yield 

' L'^ — kr' 


Q = 9 + arctan 


Lvr 


(4.28) 


Similarly, the temporal hrst integral fl3.3ip with tq = r* gives 
rv sgn{v)\L\^yk'^ — 2\E\L‘^ 


T =t + 


2\E\ 


2\E\k 


k 


(2|B|)3/= 


arctan 


/ k-2\E\ 

(tpp 


— sgn(n) arctan 


\/k^ -2 \E\lA 

VW\L ) 


(4.29) 


From the hrst integral fl4.28p . the curves describing trajectories with F^ > 0 are algebraically 
given by 


r{9) 


F2 

k + ^k‘^-2\E\L'^ sin (0 - 9)' 


(4.30) 


These curves are ellipses fl4.1ip . with a diherent angular parameterization compared to the 
turning point case. 

On each elliptic trajectory, both hrst integrals 0 and T are piecewise functions which are 
continuous at the two inertial points r = r* but have a jump discontinuity at the two apsis 
points r = rmin and r = rmax- The jump in 0 is equal to 


A0 = ^A9 = 71 (modulo 27r) (4-31) 

which is the angular separation between the two apsis points. Similarly, the jump in T is 
half of the period of a single closed orbit, 

AT = 2 At = 

The physical and geometrical meaning of these jumps is directly related to the division 
of an elliptic trajectory into two rehection-symmetric pieces around the radial line (semi¬ 
major axis) through the two apsis points. Note that the two inertial points r = r* lie on a 
perpendicular line (semi-minor axis) with respect to the apsis points. Therefore, when the 
hrst integrals 0 and T are evaluated at any point {r{t),9{t)) on the piece with sgn(u) ^ 0, 
we obtain (r(T), 9{T)) = (r*, 0) which yields the angle of the radial line that intersects the 
inertial point on that piece of the trajectory, and the time at which this point is reached. 

Next, for E = 0, the angular hrst integral fl3.30l) with tq = r* is again given by the 
expression fl4.28p . while the temporal hrst integral fl3.3ip with tq = r* yields 

r|u|(F^ -|- kr) — 2|F|^\ 

SF ) ■ 


77k 


VW\ 


(4.32) 


T = t — sgn(u) 
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(4.33) 





































(4.34) 


The resulting curves that describe trajectories with E = Q are algebraically given by 

r2 

r(0) = -. 

^ ’ fc(l + sin(0-0)) 

These curves are parabolas (I4.11I1 . with an angular parameterization that is shifted by 7r/2 
compared to the turning point case. The two inertial points on each parabola are related by 
reflection symmetry through the radial line that connects the periapsis point r = to the 
origin r = 0. 

On each parabolic trajectory, both first integrals 0 and T are piecewise functions which 
are continuous at the two inertial points r = r* but have a jump discontinuity at the 
periapsis point r = rmin, corresponding to the division of the trajectory into two reflection- 
symmetric pieces. The jumps in 0 and T are given by the same respective values fl4.31l) 
and fl4.32l) as in the elliptic case. Therefore, when the hrst integrals 0 and T are evaluated 
at any point {r(t),9(t)) on the piece of a parabolic trajectory with sgn(u) ^ 0, we obtain 
{r{T),6{T)) = (r*,0) which yields the angle of the radial line that intersects the inertial 
point on that piece of the trajectory, and the time at which this point is reached. 

Last, for > 0, the same expression fl4.28p is obtained for the angular hrst integral fl3.3UI) 
with ro = r*, while the expression for the temporal hrst integral fl3.31l) with ro = r* is given 
by 


^ rv sgn{v)\L\y/k‘^ -|- 2 EL‘^ 


k 


arctanh 




vr 


\ 2Er + k 


sgn(.)arcta„h 


(2E)3/2 

Trajectories with E > Q are algebraically described by the curves 

r2 

r{e) = 


k + sin(0 - 9) 


(4.35) 


(4.36) 


which are hyperbolas (I4.23p . with an angular parameterization that is shifted by 7r/2 com¬ 
pared to the turning point case. The two inertial points on each hyperbola are related by 
rehection symmetry through the radial line that connects the periapsis point r = r^in to the 
origin r = 0. 

The physical and geometrical meaning of the hrst integrals 0 and T is similar to the 
parabolic case. Evaluation of 0 and T at any point {r{t),9{t)) on the piece of a hyperbolic 
trajectory with sgn(u) ^ 0 yields the angle of the radial line that intersects the inertial point 
on that piece of the trajectory, and the time at which this point is reached. At the periapsis 
point r = r*, 0 jumps by the value fl4.3ip . while the jump in T is given by 


AT 


k 


arctanh 


f \ 

\Vk^ + 2EL^) ■ 


(4.37) 


4.2. Inverse square force with cubic corrections. Consider an inverse-cube correction 
to the inverse-square central force fl4.1|) . 

E{r) = —kr~‘^ — nr 
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which has the potential 


Vir) = ( 4 - 39 ) 

Both terms in this force will be attractive if A; > 0 and k > 0, in which case the potential 
has only one equilibrium point, Teq, = oo. Since U{req.) = 0, the effective potential is given 
by 

= (4.40) 

zr^ V 

Then 


E = — K,)r — kr ^ (4-41) 

is the energy first integral fl3.26p . All trajectories can be classified qualitatively from the 
equation E = t/eff.(r). 

If K > the correction term will dominate the angular momentum term. In this case, the 
effective potential will have no extrema, and consequently there are no bounded trajectories, 
while all unbounded trajectories pass through the origin r = 0. 

Thus, hereafter we will take 

0 < K < (4.42) 


In this case, the properties of the effective potential are similar to the Kepler case. There is 
one extremum, which is a negative minimum, t/™ = — As a consequence, 
all trajectories with > 0 are unbounded, while all trajectories with 0 > are 

bounded, and bounded trajectories with E = are circular. The main difference com¬ 
pared to Kepler case will be that here the unbounded trajectories are open and precessing. 
This is called the Newtonian revolving orbit problem. 

We now evaluate the angular and temporal first integrals fl3.30|) and fl3.3ip . Note these 
first integrals exist only for non-circular trajectories. We will need the relation 


n| = \/2E -|- (k — + 2kr~^. 


(4.43) 


4.2.1. Turning points for the Newtonian revolving orbit problem. The turning points r = r^ 
of the effective potential fl4.40|) are given by replacing with — k in the expressions for 
the turning points in the Kepler case. Likewise, the minimum energy becomes 


2(L2 - «) ^ °- 


(4.44) 


Only the trajectories with E > i^min need to be considered, since the trajectories with 
E = Emin are circular (in which case the first integrals 0 and T do not exist). 

For 0 > E > Emin, there are two turning points 


= 


k± ^yk^ -2\E\{L‘^ - k) 


2\E\ 


(4.45) 


The angular first integral fl3.30p with tq = can be straightforwardly evaluated to yield 


0 ± =0 + 




arctan 


Lf — K — kr 


rv 




K 


± 


sgn(n)E TT 

VL^-n2 


(4.46) 


20 























after use of the relation fl4.43l) . Similarly, the temporal hrst integral fl3.3ip with tq = 
yields 

^ rv k / I \/2\Er — k\ tt 

n = t - — + f arctan „ =F jaga(a)), (4,47) 


2\E\ (2|E|)3/2 


Y ^^2\E\rv 

From the hrst integral (14.461) . trajectories with 0 > E > i^min are algebraically described by 
curves 


r{9) = 


K 


(4.48) 


k T ^F-2|E|(L2-k) cos(v/ 1 - K/h2 (0± - e)) 

which are bounded and have r^ax = and rmin = r*_. The angular domain of these curves 
is given by 

+ Q±<e < + 0± (4.49) 

where, as shown by Proposition [31 

ttL 

> TT (4.50) 


^A6 = sgn(n)(0+ - 0_) = 


\/L^ — K 

is the angular separation between the two turning points. Since A6 is greater than 27r, the 
curves are composed of pieces, such that the angular separation 


Ae = 


between two successive 

periapsis points r = r*_ or two successive apoapsis points r = r*_ on adjacent pieces is given 
by 

(4.51) 




K 


Moreover, by Proposition [2l adjacent pieces are rehection-symmetric around a radial line 
connecting the origin r = 0 to either of the points r = where the pieces join. If A6'/(27r) 
is a rational number, then the curve is closed, consisting of a hnite number of pieces. If 
instead A6/{2n) is an irrational number, then the curve is open, consisting of an inhnite 
number of pieces. In either case, the curve describes a precessing bounded trajectory. 

On each bounded trajectory, both hrst integrals 0± and T± are piecewise functions. Their 
physical and geometrical meaning depends on the choice of the radial value tq = r_|_ or 
ro = r_. If Vq = r+ then (r(T+), 6*(T+)) = (r^ax, ©+) is an apoapsis point on the trajectory, 
or instead if rg = r_ then (r(T_), 6*(T_)) = (rmin, 0_) is a periapsis point. In both cases, 
Q± and T± are continuous at the apsis point r = but have a jump discontinuity at 
the opposite apsis r = r*:p [2]. The jump in 0± is equal to the angular separation fl4.51l) 
between two successive periapsis points or two successive apoapsis points on the trajectory. 
The jump in T± is the corresponding time interval between these points. 

At = (4.52) 

yw 

which represents the period of a single open orbit. Both of the temporal hrst integrals T± 
thus yield a sequence of times such that the time interval fl3.39p is the period At. 

Therefore, when the angular hrst integral 0_ and temporal hrst integral T_ are evaluated 
at any point {r(t),9(t)) on a piece of the bounded trajectory fl4.48p . we obtain the angle 
of the radial line that intersects the periapsis point on that piece, and the time at which 
this point is reached. Likewise the angular hrst integral 0+ and temporal hrst integral T+ 
evaluated at any point (r(t), 9(t)) on a piece of the bounded trajectory yield the angle of the 
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radial line that intersects the apoapsis point on that piece, and the time at which this point 
is reached. 

Next, for E = Q, there is a single tnrning point, 

— K 


r* = 


2k 


The angnlar hrst integral fl3.30p with tq = r* is straightforward to evalnate, 

L f — K — kr\ sgn{v)L vr 

EE^ — K, 2 


0 = 0 + 




arctan 


K 


rv 




K 


Similarly, evaluation of the temporal first integral fl3.3ip with tq = r* yields 

rv{L?‘ — K — kr) 


T = t 


3P 


(4.53) 


(4.54) 


(4.55) 


From the hrst integral (I4.54p . trajectories with E = 0 are algebraically described by curves 

E^ — K 


r{6) = 


k{l + cos(-^l — k/L"^ (0± — 9))) 


on the angular domain 


ttL 




0 < 0 < +0 + 


ttL 


K El? — K 

These curves (I4.56p are unbounded and have = r*. By Proposition |2l each curve is 
rehection-symmetric around the radial line connecting the origin r = 0 to the periapsis point 
r = r*, where the two pieces differ by sgn(n) ^ 0 (with sgn(n) = 0 at the turning point). 

On each unbounded trajectory, the hrst integrals 0 and T are piecewise functions such that 
(r(T),0(T)) = (r^nin, 0) is the periapsis point. Both 0 and T are continuous at this point. 
Thus, when these hrst integrals are evaluated at any point (r(t),0(t)) on the trajectory, we 
obtain the angle of the radial line that intersects the periapsis point and the time at which 
this point is reached on the trajectory. This is qualitatively the same as the Kepler case. 
Last, for > 0, there is again a single turning point, 

Ek^ + 2E{L^ -K)-k 


(4.56) 


(4.57) 


r* = 


2E 


(4.58) 


Evaluation of the angular hrst integral fl3.30p with tq = r* straightforwardly yields the same 
expression fl4.54p as in the case E = 0. The temporal hrst integral fl3.3ip with rg = + gives 


T = t- 


1 


rv + 


k 


E^ 


arctanh 




vr 


2Er + k 


(4.59) 


From the hrst integral fl4.54p combined with the relation fl4.5p . trajectories with E > 0 are 
algebraically described by curves 


r{6) = 

on the angular domain 

71 


K 


k + — 2\E\{L? — k) cos(^1 — njl? (0± — 0)) 


-arctan [ E2EEE^ — + 0 < 0 < —h arctan [ E2EEE^ — n/k? + 0. (4.61) 


(4.60) 
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The properties of these curves as well as the angular first integral 0 and temporal first 
integral T are the same as in the case E = 


4.2.2. Inertial points for the Newtonian revolving orbit problem. The effective potential 
fl4.40p has a single inertial point r = r*, which is given by replacing with — k in 
the expression fl4.26p for the inertial point in the Kepler case, 


r 


* 


— K 


(4.62) 


All of the earlier discussion of the first integrals 0 and T using tq = r* in the Kepler case 
carries over here in the same way as the discussion using turning points. In particular, 
the main qualitative difference compared to the Kepler case is that for 0 > > E mir, the 

trajectories are composed of more than two pieces, so consequently 0± is multi-valued even 
when the trajectories are closed. 


5. Symmetry formulation 


The general form of Noether’s theorem provides a one-to-one explicit correspondence 
between hrst integrals and dynamical symmetries of a Lagrangian. This correspondence 
arises from two variational identities as follows. 

We start with the polar Lagrangian fl2.4l) for the central force equations of motion fll.4p . 
using the polar variables t, r, 6 , v, uj. Now consider any vector held fl2.9l) . with the char¬ 
acteristic form fl2.11ll . given by components T(t,r,9,v,u;), f(t,r,9,v,u;), and 'if(t,r,9,v,u;). 
The action of this vector held on the Lagrangian can be expressed in terms of the equations 
of motion, through the variational derivatives (12.61) after integration by parts, yielding the 
identity 

prX(L) = + P^^ + ^, ^ = P^L, + P^L^. (5.1) 

or o9 at 

From this identity, a necessary and sufficient condition for the Lagrangian to be invariant 
modulo a total time derivative is that the components P^ and P® of the vector held have to 
satisfy the condition 


5h 

6 r 69 dt 


(5.2) 


for some function R{t, r, 9, v,u). A vector held fl2.1ip having this property is called a varia¬ 
tional symmetry. Because any total time derivative is annihilated by a variational derivative, 
the extremals of the Lagrangian are preserved by variational symmetries, and hence every 
variational symmetry is an inhnitesimal symmetry of the equations of motion. Thus, varia¬ 
tional symmetries have the characterization as inhnitesimal symmetries of the equations of 
motion (II.4p such that the invariance condition 


prX(L) = S = (5.3) 

holds for some function R{t,r,9,v,u). Note that a variational symmetry will be an inhn¬ 
itesimal point transformation if (and only if) r, have no dependence on v and 9, or 
equivalently, P^ and P® satisfy 

PEv = 0, P'.. = 0, PE = 0, P\ = 0. (5.4) 
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To connect the invariance condition fl5.3p to first integrals, consider the time derivative of 
a general hrst integral I{t, r, 6, v,uj) of the central force equations of motion fll.dp . The chain 
rule combined with the determining equation fll.51) directly yields the identity 


,-2; 


(5.5) 


dl _ _Sh _ 5L 

dt 6r 66 

Note the variational derivatives will vanish precisely for solutions of the equations of motion. 
The pair of coefficients of these derivatives are called the multiplier of the first integral, which 
we will denote 


Q'- = - 4 , = 


(5.6) 


From determining equation fll.5p . it is simple to see that there are no first integrals depending 
only on the variables t, r, 6. Consequently, any two first integrals that differ by at most a 
constant will have the same multiplier, as given by equation (I5.6p . Conversely, any multiplier 
determines a first integral to within an additive constant through inverting equation fl5.6l) 
by a line integral, 

1 = f {-Q^dv-r^Q^du) (5.7) 

Jc 

where C denotes any curve in the coordinate space (n, ca), starting at an arbitrary point and 
ending at a general point (n,a;). 

A comparison of identities fl5.5p and fl5.ip combined with the invariance condition fl5.3p 
now leads to the following general form of Noether’s theorem. 


Proposition 4. For the polar equations of motion fll.4p . variational symmetries in charac¬ 
teristic form X = and multipliers (Q^, Q^) for first integrals have a one-to-one 

correspondence given by 

r = Q'-, = Q® (5.8) 


and 

R = I + const. (5.9) 

In particular, through the relations (15.61) and (15.8|) . variational point symmetries correspond 
to first integrals that are at most quadratic in v and u. 


From this result, it is straightforward to obtain the symmetries that correspond to the 
first integrals (13.251) . (I3.26p . (I3.30p . (I3.3ip admitted by the polar equations of motion. 

The first integrals (13.251) for angular momentum and (I3.26p for energy have the respective 
multipliers 

(0", = (0, -1), (Q'-, Q’)iE) = (-V, -w) (5.10) 

which yield the infinitesimal symmetries (in characteristic form) 

X(L) = -do, X(E) = -vdr - ude. (5.11) 

From the relation (12.lip , these symmetries clearly correspond to infinitesimal point trans¬ 
formations 

X(L) = —Oq, X(£;) = dt- (5.12) 

Compared with the point transformations found earlier in Theorem [T] for a general central 
force, we see that —X(i) = Xi is the generator of a group of rotations (I2.2UI) and '^{e) = X 2 
is the generator of a group of time-translations (I2.2ip . 
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(5.13) 


The multiplier of the first integral fl3.30p for the angular quantity 0 is given by 

Q^){e) = {—vQe, —©l — ojQe)- 
This yields the inhnitesimal symmetry (in characteristic form) 

X(©) = —V^Edr — (©L +i^©e)56»- (5-14) 

Since © depends on L and E which themselves depend on the variables v and cj, the symmetry 
X( 0 ) is not an inhnitesimal point transformation (namely, its components do not satisfy 
condition fl5.4p b However, this symmetry can be converted into a simpler form by putting 
r = ©E in the relations fl2.1ip and fl2.9p . giving the equivalent dynamical symmetry 

X(e) = QEdt - QLde. (5.15) 

The hrst integral fl3.3ip for the temporal quantity T has the analogous multiplier 

(Q^ Q^)(t) = {-vTe, -Tl - ujTe). (5.16) 

This yields the inhnitesimal symmetry (in characteristic form) 

X(T) = —vTEdr — {Tl + ujTE)dd (5-17) 

which, similarly to X( 0 ), is not an inhnitesimal point transformation (since its components 
do not satisfy condition fl5.4p b A simpler equivalent symmetry is obtained by putting t = Te 
in the relations fl2.1ip and fl2.9p . yielding the dynamical symmetry 

X(T) = TEdt — Tide- (5.18) 

Proposition 5. Through Noether’s theorem for the polar equations of motion fll.4p . the first 
integrals for angular momentum fl3.25p and energy fl3.26p correspond to the infinitesimal 
point symmetries fl5.12p . while the first integrals given by the angular quantity fl3.30p and the 
temporal quantity fl3.3ip correspond respectively to the infinitesimal dynamical symmetries 
flHT^ and HATSD . 

The inhnitesimal dynamical symmetries given by X( 0 ) and ^{t) each generate a group 
of transformations acting on the dynamical variables {r{t),6{t),v{t),uj{t)) for any solution 
(r(f), 0{t)) of the polar equations of motion fll.4p . We will see later that there is a simple way 
to derive the transformations by hrst looking at how X( 0 ) and X((r) act on the hrst integrals 
L and E. This will also lead to a simple way to hnd the structure of the four-dimensional 
group of transformations generated by all of the inhnitesimal symmetries X(i), X(^), X( 0 ), 
and X((r). 

5.1. Transformation of first integrals under dynamical symmetries. We now work 
out how the inhnitesimal symmetries X(i), X(e), X( 0 ), and ^{t) act on all of the hrst 
integrals L, E, ©, and T. Since every hrst integral is a function of the variables t,r,6,v,u}, 
the symmetries need to be prolonged to the coordinate space {t, r, 9, v, u). This is most easily 
carried out through the prolongation relations fl2.13p and fl2.12|) which use the equivalent 
symmetries X(i), X(^), X( 0 ), and For a general dynamical symmetry given by a 

generator fl2.9p . its action on a hrst integral I{t,r,9,v,u) is thus given by 

P^X(^)Loin. = (P’^X(/) + 1^^^^ = prX|,oi„.(/) 
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which involves only the prolongation of the characteristic generator fl 2 . 1 ip evaluated on 
solutions of the polar equations of motion fll.4p . 

In the case of the infinitesimal point symmetries ^{l) and ^{e), the prolongation of their 
characteristic generators fIS.lip evaluated on the equations of motion fll.dp is straightfor¬ 
wardly given by 

P’^^wlsoin. = ~ + (2w'i;/r)a^. (5.20) 

A direct calculation (using Maple) then yields 

P’'^wlsoln.(^) = P’'^(®)lsoln.(^) = 0’ = O' 

This result means that L and E can be treated as constants under the action of the symme¬ 
tries X(/,) and X(E). Hence the action of these symmetries on 0 and T becomes simple to 
calculate, yielding the result 

P^^(0Loln.(^) = P^^(^)lsoln.(®) = 0’ P^^(0lsoln.(®) = I soln. (^) = “I' (5-22) 

The prolongation of the characteristic generators fl5.14p and fl5.17p of the infinitesimal 
dynamical symmetries X( 0 ) and X(t) are more involved to calculate: 

P^^(0)lsoln. =^(0) + 

_.2 / ^<9 

V dt ^ dt ^ ^v/r)~Ej a;, 

=^T) + 0 , 

\ ^ (5-24) 

, [ dTi I dTE I / \T \ f) 

One simplification is that, since L and E are constants of motion, the time-derivative terms 
can be calculated just by using the chain rule expression fl2.14p . This leads to the simple 
result (using Maple) 

= 0, prX(0)l^^j^ (L) = -prX(r)l^^i^ (^) = 1- (5-25) 

From these expressions, it is straightforward to calculate the action of the symmetries X( 0 ) 
and X(T) on 0 and T. The final result (using Maple) is again simple, 

P’'^(0)lsoln.(©) = P’'^T)lsoln.(©) = 0’ P^^l©) Loin. (^) = I soln. (^) = 1' (5-26) 

By applying equation fl5.19p to the previous expressions, we now have the following useful 
result. 

Theorem 3. For the first integrals L, E, 0, and T of the polar equations of motion fll.4p . 
the action of the infinitesimal point symmetries X(i), X(e) consists of 

P>'X(L)l„,..(i) = Pi-XidLiJS) = prXioL,, (T) = 0, prX,„|^,je) = -1, (5.27) 

P'^X(E)1„,„ (i) = (B) = prX,E,l^__,__ (0) = 0, prX,E) 1„,„ (T) = 1, (5.28) 
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and the action of the infinitesimal dynamical symmetries X(e) and ^{t) consists of 


P^X(e)|soln.(^) = P’'X(e)lsoln.(©) = P^X(0) Loin. (^) = 0’ P^^^Ce) Loin. (^) = (5-29) 

P’'^mLoln.(^) = P^^mLoln.(®) = P^^mLoln.(^) = 0’ P^^mLoln.(^) = “I' (^'SO) 

Hence the four first integrals are canonical coordinates for these four infinitesimal symme¬ 
tries. 

We will next use this result to give a simple derivation of the group of transformations 
generated by each of the inhnitesimal dynamical symmetries X( 0 ) and ^{t)- 

The components of these two symmetries are functions only of r, L and E, where L and 
E are the hrst integrals fl3.25p and fl3.26p . Hence, the group of transformations can be found 
by integrating a system of equations involving only the basic variables t, r, 6, along with L 
and E regarded as auxiliary variables. From the expression flS.lSp for X( 0 ), this system is 
given by 


= OEirie), L(£), E{i)), ^ = eUr(e), He), E{e)) 


de 


dr{e) 


= 0 , 


dL{e) 


de 

= 1 , 


dE{e) 


= 0 


de ' de ’ de 
where e is the group parameter. Integration of equation fl5.32p easily gives 

r(e) = r, L(e) = L + e, E{e) = E. 

Equation fl5.3ip then becomes 




dt{e) 


de 

= & e {' c , T + e, E) = 


de 

(L + e)0i;,(r, L + e, E) 
2{E-U{r)) 

dQo{r, Le, E) L + e 


(5.31) 

(5.32) 

(5.33) 

(5.34) 

(5.35) 


de 2{E-U{r)y 

which can be straightforwardly integrated. This yields 

Olyd) = 6 -\- 0 o (?"5 L, E) — 0 o (^5 T + e, E) = 0 — 0 o ('^) T + e, E) 
L{Qo{r,L,E)-Qo{r,L + e,E)) 


t(e) — t + 


(5.36) 


2{E-U{r)) 


which is the transformation group generated by X( 0 ). Similar steps applied to the expression 
fl5.18p for X('r) leads to the transformation group 


with 


t{e) = t + To(r, L, E) - To(r, L, E - e) = T - To(r, L, E 
^ I r^(To(r,L,E)-To(r,L,E-e)) 


r(e)=r, L{e) = L E{e) = E - e. 
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Recall that the respective transformation groups generated by X(i) and X(e) are simply 


e{e) = 9-e (5.39) 

t{e) = t + e (5.40) 


with 

r(e) = r, L{e) = L E{e) = E. (5.41) 

Finally, we can use the same approach to obtain the Lie algebra generated by the inhnites- 
imal dynamical symmetries X( 0 ) and ^{t), together with the inhnitesimal point symmetries 
X(i) and X(£;). Through Theorem [31 we extend these symmetries to act on L and E as 
auxiliary variables: 

Xll] = -de, ^1% = du (5.42) 

X^^* = SEdt - OLdg + Ol, X^^* = TEdt - ndg - dE. (5.43) 

This extension can be viewed as a prolongation that is restricted to the variables (t, r, 0, L, E) 
appearing in the hrst integrals. 

The extended generators fl5.42p - fl5.43p have the form of inhnitesimal point transformations 
on the space {t,r,6, L, E). A straightforward calculation (using Maple) shows that the 
commutators of these generators vanish. Thus, we have the following new result. 


Theorem 4. For the polar equations of motion fll.4p . the the infinitesimal dynamical sym¬ 
metries X( 0 ) and X(T), and infinitesimal point symmetries '^{l) and together generate 

a four-dimensional abelian group of transformations fl5.36p . fl5.37p . fl5.39p . fl5.40p . When 
these transformations are extended in the form fl5.33p . fl5.38p . and fl5.4ip . then the extended 
group acts as commuting point transformations on {t, r, 6, L, E) under which r is invariant. 


5.2. Method of extended point symmetries for finding first integrals and hidden 
dynamical symmetries. The results in Theorem 0] and Theorem [3] have the important 
consequence that we can now formulate an algorithmic symmetry method to obtain all of 
the hrst integrals along with their underlying hidden symmetry group by using only point 
symmetries combined with Noether’s theorem. The method has four steps. 

First, we hnd the variational point symmetries of the polar equations of motion fll.4p . 
and through Noether’s theorem, we obtain the corresponding hrst integrals. This will give 
X = de and X = dj, yielding L and E. For later use, we also identify the joint invariant of 
these point symmetries in the space {t,r,6), which is r. 

Second, we re-write the polar equations of motion in the form of an equivalent hrst-order 
system using only the basic dynamical variables r, 6, plus the two hrst integrals, L, E (after 
elimination of v and uj in terms of L and E). This system consists of the dynamical equations 


dL 

dE ^ 

dO 2 t 

dr 

^ = 0 , 

— = 0 , 

— = r-^L, 


dt 

dt 

dt 

dt 


±^2E -r-‘^L‘^ -2U{r). 


(5.44) 


Third, we hnd point symmetries of this hrst-order dynamical system under the restrictions 
that they commute with the variational point symmetries of the original equations of motion 
and that they preserve the joint invariant of those symmetries. The form of the symmetry 
generators being sought is thus 


X = pfir, L, E)dt + v\r, L, E)dg + r/^(r, L, E)dL + p^{r, L, E)dE. 
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This will give the four symmetries fl5.42p . fl5.43p . up to expressions which are arbitrary 
functions of L and E. 

Last, we find the canonical coordinates of these point symmetries in the space (t, r, 6 *, L, E). 
This will reproduce the first integrals L and E as well as yield the two additional first integrals 
0 and T. 

It is important to emphasize that the determining equations for point symmetries with 
the special form 05.451) will be a linear system of four equations in four unknowns, whose 
solution will not involve the hrst integrals 0 and T. In particular, the system simply consists 
of 

drlTj^ = 0, drf]^ = 0, (5.46) 

, ^ r-^Lr]^ - 7]^ Q ^ r-\2{E -U{r))'q^ - Lr]^) 

" ' {2{E - U{r)) - r-‘^L^fn ’ " ' (2{E - U{r)) - ^ > 

This is in contrast to the case of general point symmetries, where the determining equations 
will be an underdetermined system whose solution will necessarily contain expressions which 
are arbitrary functions of all of the first integrals. 


6. First integrals in n dimensions 


We return to the central force equations of motion in the n-dimensional form fl 2 .ip given 
in terms of the position vector f and velocity vector v in M”. Hereafter, we will put m = 1 
without loss of generality (via rescaling the physical units of the dynamical variables). 

We recall that all solutions {f{t),v{t)) of the equations of motion he in a time-independent 
plane spanned by these two vectors. Let f be the unit vector along r and 9 be the unit vector 
orthogonal to Fin the plane of motion, such that {f, 6} is a right-handed orthonormal basis 
in this plane C M”. Then the vectors f and v are related to the polar variables (r, 6, v, u) 
by 

f=rf, V = vf + ru6 ( 6 - 1 ) 

and 

Cl = cos( 6 *)f — sin(6')0, 62 = sin( 6 *)f -f cos{9)9 ( 6 . 2 ) 

where { 61 , 62 } is a hxed (time-independent) orthonormal basis aligned with the respective 
directions given by 6 ^ = 0 and 0 = 7 r /2 in the plane of motion. Thus we have 


r = \r\ = r ■ f, v = r ■ v, 00 = r ^6 ■ v 


and 


cos( 6 *) = f ■ ei = 9 ■ 62 , sin( 6 *) = f ■ 62 = —9 ■ ei. 

From the equations of motion fl2.ip for we easily derive 

df ^ d9 

— = uj9, — = —ujr 
dt dt 

which implies that the antisymmetric product of f and 9 is time-independent 

L(f A«) = 0. 


(6.3) 

(6.4) 

(6.5) 


( 6 . 6 ) 


Thus 


h = f A 9 
29 


(6.7) 










( 6 . 8 ) 


is an antisymmetric tensor (bi-vector) in M”, satisfying 

dt . 

—;— ~ 0, L A L = 0, L = Cl A 02 
dt 

and 

f-t = e, §-t = -f, L ■ L = 2. (6.9) 

There is a one-to-one correspondence between this tensor and the plane of motion. It fol¬ 
lows that L is a hrst integral of the n-dimensional equations of motion and determines the 
orientation of the plane of motion. Note this hrst integral does not appear among the polar 
hrst integrals since it has no dynamical content for the motion within that plane. 

Using these preliminaries, we now will express the polar hrst integrals fl3.25p . fl3.26p . fl3.30p . 
fl3.3ip in a geometric n-dimensional form. 

We start with the energy fl3.26l) . We re-write this hrst integral by noting 

|np = -I-(6.10) 

so thus 

E = + U{\f\) (6.11) 

is directly in an n-dimensional form. 

Next we consider the angular momentum fl3.25p . This hrst integral can be written as 

L=\f\e-v (6.12) 


by using the vector 9. A more physically and geometrically natural formulation is given by 
combining L and L into the hrst integral 


L = LL = Lr A9. 

From relation fl6.1l) . we can express 

9 = L~^{\r\v — {v ■ r^f) 

in terms of f and v, yielding 

L = r A F. 

This hrst integral is an antisymmetric tensor (bi-vector) in M” having the properties 


dL 

dt 


= 0, |Lp = 2|Lp, LAL = 0. 


We next turn to the angular quantity fl3.30l) . which can be expressed as 

0 = 0 + $ 

in terms of the integral expression 

sgn(r>) 


$ = -L 


'ro 


^2{E + f/(req.) - U{r))r^ - L‘^r^ 


dr. 


(6.13) 

(6.14) 

(6.15) 

(6.16) 

(6.17) 

(6.18) 


Since 0 is an angle in the plane of motion in it can be geometrically formulated as a 
unit vector 

0 = cos(0)ei -|- sin( 0 )e 2 . 
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(6.19) 



















( 6 . 20 ) 


This vector is a first integral and has the properties 

^=0, |©| = 1, ©AL = 0. 

at 

We can write it in terms of the unit vectors f and 6, or equivalently the position and velocity 
vectors r and v, by substituting the expressions (16.17^ and (I6.2p . followed by using the 
relations fl6.3p and fl6.14l) . which gives 

-'1/ T ■ V 

© = cos(<h)r + sin($)6' = — cos(<h)-sin(<h) r + — sin(<h)-y (6-21) 

|r| V L / L 

where the integral expression $ involves L, E, |f^, and also sgn(n). Note the dependence on 
sgn(n) drops out of cos(<h) = cos|<h| (since it is an even function). It is useful to re-write 
— sin(<f)) = sgn(n)sgn(L) sin |<f)| so that sgn(n) appears as an overall factor. After re-writing 
these terms in the expression cos(<h)r + sin(<f))0, and using the relation fl6.9p . we obtain 


© 


|r| M cos(0o)r — \/2sgn('y • r)|L| ^sin(0o)f^'L 


( 6 . 22 ) 


which has a geometrically simple form, where 


00 = 



sjA{E + U{r,^)-U{r))i 


--1 


12.^2 


dr 


(6.23) 


is a time-dependent angle. Here req, is an equilibrium point fl3.22p of the central force 
potential U{\f\), and tq is a turning point fl3.29p or an inertial point fl3.28p of the effective 
potential 

(7.,.(|f|) = |L|V(2|r1)" + I7(|f|) - I7(r.,.). (6.24) 

We remark that another vector first integral can be formed from 0 by rotating © by the 
angle 7r/2 in the plane of motion. This yields 


©■‘- = — sin(0)ei -I- cos(0)e2 = — sin($)r -|- cos(<h)0 (6.25) 


which we can write equivalently in the geometrical form 

©■'■ = © ■ L = (|?^L)“^ ^cos(0o)h’- L -I- )^sgn(-y • r)|L| sin(0o)f^j • (6.26) 

It has the properties 


d&^ 

dt 


0 , 




1, 0^AL = O, ©^-0 = 0. 


(6.27) 


Finally, we consider the temporal quantity (13.3ip . This hrst integral has the n-dimensional 
form 

T = t — \/2sgn('i; • r)ro (6.28) 


with 


"Id 


A = 


'ro 


A(B+(7(r„.)-C/(r)) 


- 1 

|LPr“^ dr 


(6.29) 


which we obtain through the relations (16.3p . Note tq is a time-dependent geometrical ex¬ 
pression, where rgq. and tq are the same radial points used in the hrst integral ©. 
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In summary, the four polar first integrals fl3.25p . fl3.26p . fl3.30p . fl3.3ip give rise to corre¬ 
sponding n-dimensional first integrals fId.lSp . flG.lip . fl6.22p . fl6.28p which each have a geo¬ 
metrical form given in terms of the position and velocity vectors r and v. Note that the 
energy 03.261) and the temporal quantity 03.3ip are scalar hrst integrals, whereas the angular 
momentum 06.15p is a tensor hrst integral and the angular quantity 06.22p is a vector hrst 
integral. We will show later that this vector hrst integral can be used to dehne a general 
Laplace-Runge-Lenz vector. 

The total number of time-independent quantities dehned by all of these n-dimensional 
hrst integrals is n(n — l)/2 from the components of L, and n from the components of © 
(as dehned with respect to any hxed orthonormal basis of R"'), plus 2 given by E and T. 
However, not all of these quantities are independent, due to the algebraic properties LAL = 0, 
0 AL = 0, |©| = 1. In particular, L has only 2n —3 independent components, corresponding 
to the orientation of the 2-dimensional plane of motion in R" (which accounts for 2n — 4 
components) and the magnitude of the angular momentum in this plane (which accounts for 
a single component); © has n — 1 independent components, corresponding to a 2-dimensional 
cone in R”' (which accounts for n — 2 components) and an angle of a unit vector in this cone 
(which accounts for a single component), where the cone is tangent to the plane of motion. 
Consequently, together L and © determine a total of 2?7, — 2 independent quantities. Since 
there are 2 independent quantities given by E and T, this yields 2n independent quantities 
altogether, which is precisely the number of functionally independent hrst integrals allowed 
for a second-order system with n independent dynamical variables. 

Thus we have the following result. 

Theorem 5. For the equations of motion fll.3p of general central force dynamics in n > 1 
dimensions: 

(1) L and E are well-defined first integrals for all solutions r{t). 

(2) © and T are well-defined first integrals for all non-circular solutions r{t). 

(3) Altogether, L, E, © yield 2n — 1 functionally independent constants of motion, while T 
is a first integral which is not a constant of motion. 

(4) Every first integral is a function of L, E, ©, T, and every constant of motion is a 
function o/L, E, ©. 

We emphasize that each of L, E, ©, T can be directly verihed to obey the dehning equation 
fll.2p for hrst integrals of the central force equations of motion fll.3p . We also emphasize that 
the integral expressions fl6.23p and fl6.29p appearing respectively in L and T are well-dehned 
because we have specihed the endpoint value tq directly in terms of the ehective potential 
fl6.24p . If the endpoint were omitted or allowed to be arbitrary, then the resulting hrst 
integrals would be dehned only to within an arbitrary additive constant and would lose an 
important part of their physical meaning, which we will discuss next. 

6.1. General Laplace-Runge-Lenz vector. The vector hrst integral © has physical prop¬ 
erties similar to the angular hrst integral 0 discussed in Sec. 13.31 These properties depend 
on the particular choice made for the radial endpoint ro in the integral expression fl6.23p 
appearing in ©. 

There are two diherent possibilities for rg, which are tied to the shape of trajectories in 
the plane of motion for solutions r{t) of the central force equations of motion fll.3p . Recall, 
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for a given non-circular trajectory, a turning point is a radial value at which the radius |r| 
of the trajectory is a local extremum (i.e. an apsis), and an inertial point is a radial value 
at which the radial speed f • n is a local extremum. For all central force potentials U{\f\) of 
physical interest, at least one turning point or one inertial point can be assumed to exist on 
each non-circular trajectory. 

When a trajectory has no turning points, then we choose tq to be the radius of an inertial 
point on the trajectory in the plane of motion. In this case, © will be the direction vector 
of the radial line that intersects the trajectory at this point. Similarly, the first integral T 
will be the time t at which the trajectory reaches this point. 

When a trajectory instead has at least one turning point, then we can choose tq to be 
this turning point on the trajectory in the plane of motion. Then © will be the direction 
vector of the radial line that intersects the trajectory at this point, and the first integral T 
will be the time t at which the trajectory reaches this point. If the trajectory possesses a 
single turning point, this point will be a periapsis point (i.e. a local minimum of the radius 
|r|), and © is then uniquely defined. In contrast, if the trajectory possesses multiple turning 
points, then the trajectory will have apses that occur in pairs consisting of a periapsis point 
and an apoapsis point (i.e. a local maximum of the radius |r|). In this case, © will be 
single-valued only on each part of the trajectory between pairs of successive periapsis points 
or successive apoapsis points. This is caused by the sign factor sgn(F • f) appearing in ©, 
which produces a jump in the value of © when the trajectory passes through the next apsis 
point occurring after the apsis point that corresponds to the radius vq. From Proposition |3l 
this jump is equal to the angular separation fl3.38p between each pair successive periapsis 
points or successive apoapsis points on the trajectory. 

If this angular separation fl3.38l) is a rational multiple of 27r, then the shape of the trajectory 
is closed. In this case the first integral © will yield a hnite number of distinct vectors in the 
plane of motion, corresponding to the finite number of periapsis points or apoapsis points. 
When the angular separation is exactly 27r, there will be exactly two apsis points and then 
the trajectory is not be precessing in the plane of motion. Otherwise, there will be more 
than two apsis points, and then the trajectory is precessing. 

In contrast, if the angular separation fl3.38p is an irrational multiple of 27r, then the curve 
describing the shape of the trajectory is open. In this case the first integral © will yield 
an infinite number of distinct vectors in the plane of motion, corresponding to the infinite 
number of periapsis points or apoapsis points, and the trajectory is then precessing in the 
plane of motion. 

In either case, the temporal first integral T yields a periodic infinite sequence of values, 
which are the times at which the trajectory passes through successive periapsis or apoapsis 
points. 

Based on these properties, it is natural to view the vector © as defining the directional part 
of a general Laplace-Runge-Lenz vector for arbitrary central forces in n > 1 dimensions. We 
first re-write the expression fl6.2ip for © in the following way, which will make contact with 
the usual Laplace-Runge-Lenz vector expression for an inverse-square central force [ami. 

From expression fl6.23p . we have the relation cos(0o) = sin(0o) which 

yields 

\t\v * T 

sgn(F- sin(0o) =-cos(0o)- 
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(6.30) 













(6.31) 


By using expressions fl6.ip and fl6.13p . we also get the relation 


(n • r^f ■ L = \r^v ■ L + L'^f. 


These relations fl6.3ip and fl6.30p can be combined with expression fl6.22p . yielding 


^ a|^^(|h1cos(0o))^ , 2|rl2^p^(cos(0o))^ ^ 

0 = —- r A - ~ V 2 - ^ 



(6.32) 


Thus, we obtain the general Laplace-Runge-Lenz vector 



(6.33) 


where A{E, L) > 0 is an arbitrary normalization factor, and 0o is the (time-dependent) angle 
expression fl6.23p . This vector satishes 


dA T* 

— = 0, A A L = 0. 


(6.34) 


From the preceding discussion, we have the following main result. 

Theorem 6. For all non-circular solutions of the equations of motion fll.31) of general cen¬ 
tral force dynamics in n > 1 dimensions, the vector quantity (16.331) is a well-defined first 
integral A whose geometrieal properties are determined by the radial value ro > 0 appearing 
in the angular quantity 0o- 

(1) A lies along the direction of the radial vector r{T) in the plane of motion of the trajectory 
r{t), where T is the temporal first integral (I6.28p using the radial value vq. 

(2) For an unbounded trajectory r{t) with no radial turning points, T yields the times t when 
|r(t)| = ro is the radius of the two inertial points on the trajectory (where the radial speed 
v(t) ■ r is a positive or negative extremum); thereby A is double-valued. 

(3) For an unbounded trajectory r{t) with one radial turning point, T yields the time t when 
|r(t)| = ro is the radius of the periapsis point on the trajectory (where the radius |r(f)| is a 
minimum); thereby A is single-valued. 

(4) For a bounded trajectory r{t) with multiple radial turning points, T yields the times t 
when |r(t)| = ro is the radius of either the periapsis point(s) or the apoapsis point(s) on the 
trajectory (where the radius |r(t)| is, respectively, either a minimum or a maximum); there¬ 
fore, A is single-valued if the angular separation between successive periapsis and apoapsis 
points is 71, and otherwise A is multi-valued. 

An interesting variant of the general Laplace-Runge-Lenz vector A can be obtained by 
choosing the radial value ro to be the radius of an inertial point for any non-circular tra¬ 
jectory. As we will see in the next section, this variant generalizes Hamilton’s eccentricity 
vector PI which is known to be a hrst integral for inverse-square central forces. 

7. Examples of u-dimensional general Laplace-Runge-Lenz vector 

We will now examine the general Laplace-Runge-Lenz vector (16.331) and its variant for two 
important examples of central force dynamics in n dimensions, using the results from Sec. 01 
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For comparison with previous results in the literature on central force dynamics in n = 3 
dimensions, it will be useful to note some 3-dimensional identities relating the antisymmetric 
product of vectors and the cross-product of vectors: 

C -{A A B) = {C- A)B - {C ■ B)A = {A x B) x C (7.1) 

{AAB)-{CAD) = 2{A-B){C ■ 3) -2{A-D){B-C) = 2{A x B) ■ {C x 3) (7.2) 


7.1. Inverse-square force. An inverse-square force fl4.ip has three different types of non- 
radial trajectories, depending on the hrst integrals (16.151) and (Id.llh for angular momentum 
L and energy E\ elliptic trajectories, with 0 > parabolic trajectories, with E = 0] 

and hyperbolic trajectories, with E > 0. 

Parabolic and hyperbolic trajectories respectively have a single turning point fl4.16p and 
fl4.2ip . which is the periapsis on the trajectory. Using the radius of the turning point for the 
radial value tq to evaluate the angular expression fl6.18p . we have 


<h = arctan 


3 ^) 


(7.3) 


where L = |L • L is the angular momentum scalar and v = v ■ f is the radial speed. This 
yields 


cos(<F) = 


kr 


rV2ElX+¥ 

Hence the vector hrst integral fl6.32l) is given by 
1 


= cos(0o), sin($) = 


—Lv 


V2EIA+W 


(7.4) 


0 = 


r^V2EL^ + P 


((L^ — kr)f— {v ■ r)r ■ L) = 


{{k/\r\)r + V ■E) (7.5) 


^U|L|2 + A:2 

in terms of the energy E and the angular momentum L. The expression fl6.33p for the general 
Laplace-Runge-Lenz vector thereby yields 


A* = - {{k/\r\)r + v ■ L) 


(7.6) 


where we have chosen the normalization factor to be A{E, L) = E\L\^ + k?- (Note the 

subscript * indicates the use of a turning point for tq in dehning the vector.) This vector 
fl7.6p lies along the radial line connecting the origin to the periapsis point on the trajectory. 

Elliptic trajectories have two turning points fl4.8p . which are the periapsis and apoapsis 
on the trajectory. Evaluating the angular expression fl6.18p by using the radius of the 
turning points for the radial value tq, we obtain 


= arctan 


— kr\ TT 


This yields 


cos(<I)±) = T- 


E — kr 


= cos(0o), sin(<h±) = ±- 


Lv 


rV2EUT¥ V2EL2 + fc2 

with 2EE -|- > 0 for E > Emm- Hence the vector hrst integral (16.321) is given by 

Tl /.-o , .7 ±1 


(7.7) 


(7.8) 


0+ = 


rW2EE + fc2 


((L^ — kr)f— {y ■ r)r ■ L) = 


^/W\ 


fc2 


{{k/\r\)f+v-'L). (7.9) 
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Then the expression fl6.33l) for the general Laplace-Runge-Lenz vector gives 


A^± = ± {{k/\f\)f+ V ■ L) 


(7.10) 


where we have again chosen A{E, L) = as the normalization factor. This 

vector (17.101) lies along the semi-major axis of the trajectory, such that y4*_ points toward 
the periapsis point and ^4*+ points toward the apoapsis point. 

Notice the relationship A^_ = A^. This means that the expression fl7.6l) provides a general 
Laplace-Runge-Lenz vector that points toward the periapsis point for all three types of 
traiectories. In n = 3 dimensions, this vector is exactly the usual Laplace-Runge-Lenz 
vector mm [ 

A^ = V X L — k\f^~^r (7-11) 

where \A\ = \/2EL'^ + k^. 

An interesting variant of the n-dimensional general Laplace-Runge-Lenz vector fl7.6p is 
given by choosing the radial value tq to be an inertial point on a trajectory. All three types 
of trajectories each have two inertial points fl4.26l) differing by sgn(n) ^ 0 on the trajectory. 
Using the radius of the inertial points for the radial value tq in evaluating the angular 
expression fIb.lSp . we have 

/ r 2 _ \ 

<h = arctan ( — - -) (7-12) 

\ Lvr J 

which yields 


cos($) = 


\L\\v\ 


V2eIaT¥ 

The vector hrst integral (I6.22p is then given by 


= cos(0o), sin(<h) = sgn(Lr>) 


kr 


rV2EL^ + k^ 


(7.13) 


® " rV2gL" + “ kr/L^)r- L) 

= ^=2==^ ((2E + fc/|r-l)f + (1 - L) 

^2{E+ k^/\L\^)\v ■E\'' 

Hence the expression fl6.33p for the general Laplace-Runge-Lenz vector gives 


(7.14) 


A* 


{{2E + k/\r\)r + (1 - 2k\r\/\L\^)v ■ L) 

v-r\ 


(7.15) 


where we have chosen A{E, L) = ^/2{E + fc^/lLp). (Note the superscript * indicates the 
use of an inertial point for tq in dehning the vector.) This vector fl7.15p is related to the 
previous vector A* by 

X ■ L = sgn(n)A* (7.16) 

which represents a rotation of sgn('i; • r)7r/2 in the plane of motion. 

Thus, the variant Laplace-Runge-Lenz vector fl7.15p lies along a line that is perpendicular 
to the Laplace-Runge-Lenz vector in the plane of motion. (In the case of elliptic trajectories, 
this line is the semi-minor axis.) Most interestingly, it is double-valued since it changes sign 
around the periapsis point on a trajectory. In n = 3 dimensions, this vector is a multiple of 
Hamilton’s eccentricity vector [TB] 

e = V + k\L\~‘^\r\~^r X L = \L\~‘^L x A* = sgn(r;)|L|“^A* 
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Although the eccentricity vector e is usually defined only for elliptic trajectories, the variant 
vector A* exists for parabolic and hyperbolic trajectories as well. 


7.2. Inverse-square force with cubic corrections. For an inverse-square force with cubic 
corrections fl4.38p . there are again three different types of non-radial trajectories, depending 
on the first integrals (Ib.lbh and flb.lip for angular momentum L and energy E\ bounded pre- 
cessing elliptic-like trajectories, with Q > E > unbounded parabolic-like trajectories, 

with E = Q] and unbounded hyperbolic-like trajectories, with E > Q. 

Bounded trajectories have turning points fl4.45p given by the periapsis and apoapsis points 
on the trajectory. Using the radius of the turning points for the radial value tq to evaluate 
the angular expression fIb.lSp . we have 


where L = |L • L is the angular momentum scalar and v = v ■ f is the radial speed. This 
yields 

— K — kr 


cos($±) = =F cos 
sin(<h±) = ±sin 


VL^-k 

L 


arccos 


arcsm 


r^2E{L^ - k) + F 
vEE^ — K 

V v/2U(L2 - k) + F 


= cos((^o), 


(■ 


(7.19) 


In terms of these expressions, the vector first integral fl6.21|) is given by 

L / LF' — K — kr 


0± = =F cos 




arccos 


— sm 


r^2E{L'^ -K) + k‘^ 
/ vELF — K 


arcsm 


f ■ L 


(7.20) 


The alternative formulation fl6.32p of this first integral is most easily obtained by re-writing 
the two terms in the expression fl7.20l) by the following steps. First, we use the relations and 
(16.1 p and (16.9p to express the second term in the form 


sin(<F±)f ■ L = — sin(<F±)f -|— sin(<I)±)F • L. 
rv V 


(7.21) 


Next, we combine the coefficients of the f terms in expressions (I7.20p and (I7.2ip to get 


where 


We then obtain 


L EEv^E-EF 

cos(<F±) -I-sin(<I)±) =-cos(<F± — T) 

rv rv 


T = arctan ( . 


©+ = 


1 ( Er^E~+EF 


V \ r^ 


rv/ 


cos(<h± — T)r-I- --sin(<F±)F ■ L 
Jj 


(7.22) 

(7.23) 


(7.24) 
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which can be written ont in an explicit form with the nse of the relation fl4.43p . The 
expression fl6.33p for the general Laplace-Runge-Lenz vector thereby yields 


^ + fc/|r1) + K/\ffx 


V ■ r 


cos 




:arccos 


=F(L^ — k — k)\f\ 
r^2E{L‘^ - k) + F 


arccos 


V ■ r 


^/2{E + k/\f^) + K/\k]^ 


±L ^|r|sin 


arcsin 


(- 


V ■ r 






\/T2 - K \\f\^y2E{L^ - k) + k^ 



(7.25) 

in terms of the energy E and the scalar angnlar momentum L. (Note the subscript * 
indicates the use of a turning point for tq in defining the vector.) One possible choice for the 
normalization factor, in analogy to the Kepler case, is A(E, L) = y^2E{L'^ — k) + The 
properties of this vector fl7.25l) are, however, quite different compared to the Kepler case. 
Since unbounded trajectories are precessing, ^4*^ is multi-valued such that, on a trajectory 
f{t), y4*_ points toward the periapsis point closest to f{t) and ^4*+ points toward the apoapsis 
point closest to r(t). This means y4*± is jump discontinuous when r(t) passes through each 
periapsis and apoapsis, respectively. 

Unbounded trajectories have a single turning point (I4.53p in the parabolic-like case and 
fl4.58p in the hyperbolic-like case. In both cases, the turning point is the periapsis on the 
trajectory. Evaluating the angular expression flb.lSp by using the radius of the periapsis 
for the radial value tq, we obtain $ = Hence the vector first integral fl6.32p is given 
by © = ©_. Then A = A_ is the expression for the general Laplace-Runge-Lenz vector 
fl6.33p . This vector lies along the radial line connecting the origin to the periapsis point on 
the trajectory. 

The expression 


R* — 


^ {^/2{E + fc/|r1) + K/\r^^x 


V ■ r 


cos 


VL^ 


K 


arccos ^ 


L ^|rjsin 


L^-K-k\r\ 

\E\^2E{L^ - k) + k^ 

V ■ r\fl? 


arccos 


V ■ r 


VL^ 


arcsin 


K 


^]f\^2E{L‘^ - k) + k"^ 


^2{E + k/\k]) + 

v-L 


(7.26) 

thereby provides a general Laplace-Runge-Lenz vector that points toward the periapsis point 
closest to f(t), for all three types of trajectories. 

Similarly to the Kepler case, there is a variant of the n-dimensional general Laplace-Runge- 
Lenz vector (17.2611 . which arises from choosing the radial value tq to be an inertial point on 
a trajectory. All three types of trajectories each have two inertial points fl4.62p differing by 
sgn(n) ^ 0 on the each part trajectory. Using the radius of the inertial points for the radial 
value Tq in evaluating the angular expression (16.1811 . we have 


$ 


L 

, arctan 

^E^-k 
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E — K — kr 


rv 




K 


(7.27) 























































which yields 


cos(<l>) = cos 


x/l2 


arccos 








sin(<l>) = sgn(u) sin 


^^2E{L^ - K) + k^ 

— K — kr 


= cos(0o), 


(7.28) 


y/l? — H 


arcsm 


'2 - - Vr^2E(L2-K) + fc2 

The vector first integral fl6.22p is then given by 

L / |n|\/L2 — K 


©± = cos 


y/L'^ — K 
+ sgn(n) sin 


arccos 


V^2E(L2-k) + F. 

L'^ — K — kr 


(7.29) 




arcsm 


K 


r^2E{L‘^ - k) + fc2 


r ■ L. 


By the same steps as in Kepler case, this leads to the expression 
’f ^ {^2{E + fc/|r1) + /t/lr^x 


V ■ r 


cos 


VL^ 


arccos 






K 


+ sgn(n)L ^|r)sin 


^|7l^2E(L2-/s:) + A;2 

L 


— arccos 


V ■ r 


( U-K-k\r\ 

arcsm 


y/L‘^-K V|rl^2E(L2-K) + P 


^2{E + k/\r\) + K/\E(^ 

n-L 


(7.30) 

(Note the snperscript * indicates the nse of an inertial point for tq in defining the vector.) 
The variant vector A* is related to the previons vector A^, by a rotation throngh an angle 
sgn(n ■ r){L/\/LP‘ — k)t:/2 in the plane of motion. Hence, in contrast to the general Laplace- 
Rnnge-Lenz vector H*, this variant vector A* changes sign when f{t) passes throngh each 
periapsis on a trajectory. 


8. Concluding remarks 

The resnlts presented in Sec. [3] and Sec. |6] provide a simple direct derivation of all 2n first 
integrals for general central force dynamics in n > 1 dimensions. This derivation is based 
on solving the determining eqnations for first integrals of the eqnations of motion in a polar 
formulation and does not involve any use of symmetry. The first integrals are shown to be 
generated by an antisymmetric tensor (bi-vector) L determining the plane of motion, the 
angular momentum L in this plane, and the energy E, which are defined for all solutions r{t), 
plus an angular quantity 0 and a temporal quantity T both of which are defined only for 
non-circular solutions r{t). The quantities L, L, E comprise 2n — 2 functionally independent 
first integrals, which constitute a complete set for all circular solutions, while the quantities 
0, T together with L, L, E comprise 2n functionally independent first integrals, which 
constitute a complete set for all non-circular solutions. 

The angular quantity 0 is used to define a general Laplace-Runge-Lenz vector whose 
geometrical and physical properties are discussed in detail in Sec. O In particular, 0 is 
shown to be single-valued if a trajectory r{t) has at most one apsis (turning point) in the 
plane of motion. In this case the general Laplace-Runge-Lenz vector has a unique direction 
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aligned with the apsis on the trajectory. If instead a trajectory has multiple apses in the 
plane of motion then 0 is shown to be single-valued when it is evaluated between pairs 
of successive periapsis points or apoapsis points on the trajectory. In this case, 0 dehnes 
a unique direction only if the angular separation between any pair of successive periapsis 
points or apoapsis points on the trajectory is an integer multiple of 27r, corresponding to the 
trajectory being closed. If the angular separation is not an integer multiple of 27r, then the 
direction dehned by 0 undergoes a discontinuous jump when the next successive periapsis 
or apoapsis point is reached on the trajectory, corresponding to the trajectory being either 
open or processing. 

A variant of the general Laplace-Runge-Lenz vector is also dehned in Sec. |6l by using a 
different form of the angular quantity 0 based on the inertial points of a trajectory f{t). In 
this case, 0 is single-valued only on each piece of a trajectory where the radial speed has a 
dehnite sign, while at each apsis point on the trajectory, 0 has a discontinuous jump. This 
variant is a generalization of Hamilton’s eccentricity vector, which applies to general central 
forces. 

These properties are explicitly illustrated in Sec. [Hand Sec.[7]for an inverse-square central 
force, where bounded trajectories are closed and non-precessing, and for a inverse-cube 
corrections to an inverse-square central force, where bounded trajectories are precessing and 
either open or closed (depending on the correction parameter). 

A symmetry interpretation of the hrst integrals 0 and T is also obtained, using the La- 
grangian formulation of the central force equations of motion. In contrast to the well-known 
origin of L and E from point symmetries given by polar rotations and time translations, the 
quantities 0 and T are shown to arise from hidden dynamical (hrst-order) symmetries. The 
explicit transformations generated by these symmetries are derived by utilizing the action of 
the symmetries on the hrst integrals L, E, 0, T. These transformations are shown to form 
a four-dimensional abelian Lie group. This leads to a novel method for directly deriving the 
hrst integrals by use of extended point symmetries outlined in Sec. |5l 

All of these new results will be further developed in future work. We plan to apply the 
method of extended point symmetries to give a new symmetry derivation of the 2n hrst 
integrals for the central force equations of motion in n dimensions and also to obtain the 
explicit hidden symmetry group underlying the general n-dimensional Laplace-Runge-Lenz 
vector. 
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